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1 Introduction 



In classical optics, the concepts of intensity and phase of optical fields have a well-defined 
meaning. The oscillating real electromagnetic field associated with single mode, E — 
Acos{^), has a well-defined amplitude and phase. Apart from a constant factor, the 
squared real amplitude, A^, is the intensity of the field. In classical electrodynamics, 
contrary to quantum electrodynamics, there is no real need to use complex numbers to 
describe the field. However, it is convenient to work with exponentials rather than cosine 
and sine functions, and complex amplitudes of the field, E — Aexp{—i^), are commonly 
used to describe the field. The moduhis squared of such an amplitude is the intensity of 
the field and the argument is the phase. Both the intensity and the phase can be measured 
simultaneously in classical optics. In quantum optics, it was quite natural to associate 
the photon number operator with the intensity of the field and somehow construct the 
phase operator conjugate to the number operator. The latter task, however, turned out 
not to be easy. 

The first attempts to construct explicitly a quantum phase operator as a quantity 
conjugate to the number operator were made by Dirac [1927]. His idea was to perform 
a polar decomposition of the annihilation operator, similar to the polar decomposition 
of the complex amplitude performed for classical fields. It turned out later that such a 
decomposition suffers from serious drawbacks, and the phase operator introduced in this 
way cannot be considered as a properly defined Hcrmitian phase operator. Susskind and 
Glogower [1964] exposed the contradictions inherent in Dirac's polar decomposition and 
introduced, instead of the phase operator that appeared to be non-Hermitian, the opera- 
tors cos$sG and sin$sG corresponding to the cosine and sine of the phase. Unfortunately, 
these cos$sG and sin$sG operators, although being Hermitian, do not commute, so that 
they cannot represent a single phase angle. Historically, the idea to use cos$ and sin$ as 
Hermitian operators describing the phase, was first raised by Louisell [1963] in his short 
Letter, but he did not construct the explicit form of these operators. Carruthers and 
Nieto [1968], in their review paper, gave a detailed record of the problems encountered 
on the way to constructing of the Hermitian phase operator and discussed thoroughly 
the properties of cos$sG ^^^d sin$sG operators. From their analysis it became clear that 
it is the boundedness of the photon number spectrum which does not allow for negative 
values, and which precludes the existence of a properly defined Hermitian phase oper- 
ator in the infinite-dimensional Hilbert space. The difficulty in finding the form of the 
Hermitian phase operator led to the widespread belief that no such operator exists, al- 
though there were a number of ingenious attempts to construct a suitable operator within 
the infinite-dimensional Hilbert space (Garrison and Wong [1970], Turski [1972], Popov 
and Yarunin [1973, 1992], Paul [1974], Damaskinsky and Yarunin [1978], Galindo [1984a, 
1984b]). It was known (Newton [1980], Barnett and Pegg [1986], Luks and Pefinova 
[1991], Luks, Pefinova and Kfepelka [1992a]) that extension of the oscillator energy spec- 
trum to negative values allows for the mathematical construction of the Hermitian phase 
operator, but this solution was unsatisfactory because of its recourse to unphysical states. 

Recently, Pegg and Barnett [1988, 1989] (see also Barnett and Pegg [1989, 1990, 1992, 
1993], Pegg, Barnett and Vaccaro [1989], Barnett, Pegg and Vaccaro [1990], Pegg, Vaccaro 
and Barnett [1990], Vaccaro and Pegg, [1990a, 1990b, 1993], Vaccaro, Barnett and Pegg 
[1992]) have found a way out of the difficulties with construction of a Hermitian phase 
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operator. The key idea in the development of the Hermitian optical phase operator was 
abandonment of the conventional infinite-dimensional Hilbert space for the description of 
quantum states of a single field mode. They introduced, instead, a state space 'H*^'^) of 
formally finite dimension together with a prescription for taking the infinite-dimensional 
limit only after c-number expectation values and moments have been calculated. This 
idea reintroduced a symmetry to the photon number spectrum, which became bounded 
not only from below but also from above, and removed the main obstacle in constructing a 
Hermitian phase operator. An essential and indispensable ingredient of the Pegg-Barnett 
construction is the way the infinite-dimensional hmit is taken, which distinguishes it from 
the quantum-mechanical constructions based on finite-dimensional spaces that have been 
studied before (Lcvy-Leblond [1973, 1976, 1977], Santhanam and Tekumalla [1976], San- 
thanam [1976, 1977a, 1977b], Santhanam and Sinha [1978], Goldhirsh [1980]), but in 
which, when the limiting procedure has been applied for the phase operator, the original 
problems reappeared. The consequences of the limiting procedure in the Pegg-Barnett 
approach have been discussed by Barnett and Pegg [1992], and Gantsog, Miranowicz and 
Tanas [1992]. The proposal of the Pegg-Barnett approach has renewed interest in the 
problem of the proper description of the quantum-optical phase. 

Almost at the same time, Shapiro, Shepard and Wong [1989] (see also Shapiro and 
Shepard [1991]) used an alternative approach based on the quantum estimation theory 
and probabihty operator measures (Helstrom [1976]) to describe the phase properties 
of optical fields. This approach does not rely on the existence of the Hermitian phase 
operator but rather on the existence of the eigenstates of the Susskind-Glogower nonuni- 
tary exponential phase operator (Susskind and Glogower [1964]). The eigenstates of 
the Susskind-Glogower phase operators form a basis for the probability operator mea- 
sures. The Shapiro, Shepard and Wong [1989] idea has gained some popularity (Hall 
[1993], Hall and Fuss [1991], Schleich, Dowling and Horowicz [1991], Braunstein, Lane 
and Caves [1992], Braunstein [1992], Hradil [1992a, 1992b], Hradil and Shapiro [1992], 
Lane, Braunstein and Caves [1993], Jones [1993], Shapiro [1993]). It turned out, however, 
that it gives for physical states [i.e., states with a finite mean number of photons (finite 
energy and its higher moments)], the same results as the Pegg-Barnett approach after the 
limit transition to the infinite-dimensional space. The eigenkets of the Susskind-Glogower 
exponential phase operators can be used in a similar fashion as coherent states (eigenkets 
of the annihilation operator) to define the resolution of the system operators; e.g., the 
phase operator (Luks and Perinova [1991, 1993b], Brif and Ben-Aryeh [1994a, 1994b], 
Vaccaro and Ben-Aryeh [1995]. In this case the ordering of the phase exponentials is 
relevant, and, if the antinormal ordering is taken, the results agree with those obtained 
from the Pegg-Barnett formalism. 

Another way to describe the phase properties of the field is to use quasiprobability 
distribution functions. The idea behind this approach is relatively simple: to integrate 
the suitable quasiprobability distribution function, such as the Glauber-Sudarshan P- 
function, Wigner function, or Husimi Q-function over the "radial" variable and getting 
in this way a corresponding phase distribution, which can then be applied in calculations 
of the mean values of the phase- dependent quantities. Since the quasiprobability descrip- 
tion of the quantum state of the system can be in some cases associated with realistic 
measurements performed on the system, this approach to the phase problem has focused 
the attention of many authors (see §2.2). The phase distribution functions obtained by 
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integrating the quasidistributions arc different for different qiiasidistributions, and they 
are all different from the Pegg-Barnett phase distribution. The situation is even worse be- 
cause for some states of the field the P-function and the Wigner function take on negative 
values, and the corresponding phase distribution can also be negative. This means that 
such phase distributions must be used with some care, but in many cases this approach 
gives results describing properly the phase properties of the field. 

Noh, Fougeres and Mandel [1991, 1992a, 1992b, 1993a, 1993b, 1993c, 1993d, 1993e] 
(see also Fougeres, Monken and Mandel [1994], Barnett and Pegg [1993], Hradil [1993a, 
1995], Hradil and Bajer [1993]) presented an operational approach to the quantum phase 
problem. Their idea is to define phase in terms of what actually is, or can be, measured. 
They do not search for a phase operator which would satisfy some mathematical criteria, 
but start their considerations from the experimental schemes usual in classical phase mea- 
surements. Examining various measuring schemes, they identify certain operators. Cm 
and Sm-i corresponding to the measured cosine and sine of the phase difference between 
two fields. As a result, Noh, Fougeres and Mandel came to the conclusion that there 
is no unique phase operator, and that different measuring schemes correspond to differ- 
ent operators. Nevertheless, recent theoretical studies (Riegler and Wodkiewicz [1994], 
Englert and Wodkiewicz [1995], Englert, Wodkiewicz and Riegler [1995]) show that the 
intrinsic Hermitian phase operator associated with the Noh, Fougeres and Mandel ap- 
paratus can be found. The phase distribution measured in this experiment, under some 
conditions, is the radius-integrated Q-function (Freyberger and Schleich [1993], Leonhardt 
and Paul [1993a], Freyberger, Vogel and Schleich [1993a, 1993b], Bandilla [1993], Khan 
and ChaTidhry [1994]). The measurements of Noh, Fougeres and Mandel are important 
since until then only isolated phase measurements were available performed by Gerhardt, 
Biichler and Litfin [1974] (Gerhardt, Welling and Frohch [1973], and for theoretical anal- 
ysis see also Nieto [1977], Levy-Leblond [1977], Lynch [1987, 1990, 1993, 1995], Gerry and 
Urbanski [1990], Bandilla [1991], Tsui and Reid [1992], Matthys and Jaynes [1980], Walker 
and CarroU [1984, 1986]). 

Quite recently, another experimental technique, optical homodyne tomography, was 
invented and applied to measurements of the quantum state of the field (Smithey, Fari- 
dani, and Raymer [1993], Beck, Smithey and Raymer [1993], Beck, Smithey, Cooper 
and Raymer [1993], Smithey, Beck, Cooper and Raymer [1993], Smithey, Beck, Cooper, 
Raymer and Faridani [1993]) allowing quantum phase mean values to be calculated from 
the measured field density matrix. This technique opens new possibilities for quantum 
measurements. Overviews of various techniques for measuring phase distributions are pre- 
sented by, e.g., Leonhardt and Paul [1993b], Paul and Leonhardt [1994], Pefina, Hradil 
and Jurco [1994] and Lynch [1995]. 

Another interesting approach to the phase problem was presented by Bergou and 
Englert [1991]. They introduced the idea of phasors and phasor bases that can be used for 
studying possible candidates for the quantum phase operators. Different phasor bases lead 
to different phase operators, and according to Bergou and Englert [1991] extrapolation of 
the classical concept of phase to the quantum regime is not unique. 

Since the absolute phase of the single-mode field is not accessible for measurements, 
and it is always the difference with respect to a reference phase that we are forced to 
deal with in real measurements, it is tempting to define the phase-difference operator as a 
fundamental quantity describing the optical phase. Luis and Sanchez-Soto [1993c, 1994] 
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have defined a Hermitian phase-difference operator, which is in fact a polar decomposition 
of the Stokes operators for the two-mode field, and it is not the same as the difference of 
the two Pegg-Barnett operators. The difference between the two is most pronounced for 
weak fields. Ban [1991a, 1991b, 1991c, 1992] has introduced yet another phase operator, 
based on the two-mode description of the field. In recent years, many different aspects of 
the quantum phase problem have been studied (Barnett, Stenholm and Pegg [1989], Nath 
and Kumar [1989, 1990, 1991b], Chaichian and Ellinas [1990], Lakshmi and Swain [1990], 
Summy and Pegg [1990], Hradil [1990, 1993b], Luks and Perinova [1990, 1992, 1993a, 
1994], Vourdas [1990, 1992, 1993], Adam, Janszky and Vinogradov [1991], Cibils, Cuche, 
MarvuUe and Wreszinski [1991], Dowhng [1991], EUinas [1991a, 1991b, 1992], Gantsog and 
Tanas [1991d], Nath and Kumar [1991a], Orszag and Saavedra [1991a, 1991b], Paul [1991], 
Tanas [1991], Wilson-Gordon, Buzek and Knight [1991], Agarwal, Chaturvedi, Tara and 
Srinivasan [1992], Bandilla [1992], Janszky, Adam, Bertolotti and Sibilia [1992], Luks, 
Perinova and Krepelka [1992b, 1994], Ritze [1992], Smith, Dubin and Hennings [1992], 
Tsui and Reid [1992] Agarwal [1993], Agarwal, Scully and Walther [1993], Ban [1993], 
Bialynicki-Birula, Freyberger and Schleich [1993], Chizhov, Gantsog and Murzakhme- 
tov [1993], Chizhov and Murzakhmetov [1993], Daeubler, Miller, Risken and Schoendorff 
[1993], d'Ariano and Paris [1993, 1994], Jex and Drobny [1993], Luis and Sanchez-Soto 
[1993a], Stenholm [1993], Schieve and McGowan [1993], Tu and Gong [1993], Agarwal, 
Graf, Orszag, Scully and Walther [1994], Belavkin and Bendjaballah [1994], Vaccaro and 
Pegg [1994a, 1994c, 1994b], Bialynicka-Birula and Bialynicki-Birula [1994, 1995], Das 
[1994], Franson [1994], Gennaro, Leonardi, Lillo, Vaglica and Vetri [1994], Opatrny [1994], 
Sanchez-Soto and Luis [1994], Schaufler, Freyberger and Schleich [1994]), and the number 
of publications on the subject is growing steadily. Various conceptions of the quantum- 
optical phase have been described by Barnett and Dalton [1993] in a special issue of 
Physica Scripta devoted to "Quantum phase and phase dependent measurements", and 
in the same issue one can find very interesting historical facts, given by Nieto [1993], 
concerning the development of our knowledge on quantum phase. 

Nowadays, although the quantum phase is still a subject of some controversy, sig- 
nificant progress has been achieved in clarifying the status of the quantum-mechanical 
phase operator, describing the phase properties of optical fields in terms of various phase 
distribution functions, and measuring phase-dependent physical quantities. We can now 
risk the statement that, despite the existence of various different conceptions of phase, we 
are en route to unified view and understanding of the quantum-optical phase. 

It is not our aim in this review article to give a detailed account of different descrip- 
tions of the quantum phase showing their similarities and differences. We will not focus 
our attention on the quantum phase formalisms, which are interesting on their own right 
and deserve separate treatment. We shall instead concentrate on the description of the 
quantum properties of real-field states which are generated in various nonlinear optical 
processes. Nonlinear optical phenomena are sources of optical fields, the statistical proper- 
ties of which have been changed in a nontrivial way as a result of nonlinear transformation. 
Quantum phase properties are among those statistical properties which undergo nonlin- 
ear changes, and fields generated in different nonlinear processes have different phase 
properties. With the existing phase formalisms, the quantum-phase properties of such 
fields can be studied in a systematic way, and quantitative comparisons between differ- 
ent quantum-field states can be made. Using the Pegg-Barnett phase formalism and the 
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phase formalism based on the s-parametrized quasidistribution functions, we will study 
a number of both one- and two- mode field states from the point of view of their phase 
properties. 



2 Phase formalisms 

At the very beginning of quantum electrodynamics, Dirac [1927] raised the idea that the 
optical phase should be described by a Hermitian phase operator canonically conjugate 
to the number operator; that is, for the single-mode field the two operators should obey 
the canonical commutation relation 

[^,n] = -i. (2.1) 

This commutation relation implies directly the "traditional" Heisenberg uncertainty rela- 
tion 

((Anf)((A$f)>i, (2.2) 

which appeared to be wrong (Susskind and Glogower [1964], Carruthers and Nieto [1968]). 
Closer investigation of the commutator (2.1) showed that the matrix elements of the phase 
operator $ in a number-state basis are undefined (Louisell [1963]): 

{n-n'){n\^\n) = -i5„„.. (2.3) 

Since it was suggested that the problems in eq. (2.3) are due to the multivalued nature of 
Louisell [1963] introduced the operators cos$ and sin$, which should, as he suggested, 
satisfy the commutation relations: 

[cos$, h] — isin$, 

sin$,n = -icos$. (2.4) 

However, Louisell [1963] did not give the explicit form of the cosine and sine operators; 
thus, his idea did not help much in solving the phase problem. Moreover, it turned out 
that the problem expressed in eq. (2.3) was not due to the multivalued nature of <E>, but 
rather to the improper application of the correspondence between the Poisson bracket 
and the commutator. Susskind and Glogower [1964] returned to the original Dirac idea 
of polar decomposition of the creation and annihilation operators and introduced the 
exponential phase operators: 

oo 

e5$(i$sG) = Y.\n){n+1\, (2.5) 

n=0 

eSp(-i$sG) = [eSp(i$sG)]^- (2.6) 
Prom eqs. (2.5) and (2.6) one obtains: 

SEp(i$SG)exp(-i"^'sG) = 1, 

SEp(-i$sG)exp(i$sG) = 1-|0)(0|, (2.7) 
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which exphcitly shows the non-unitarity of the Susskind-Glogower phase operator. 

The Susskind-Glogower exponential operators (2.5) and (2.6) allow for constructing 
two Hermitian combinations corresponding to cosine and sine of the phase. However, the 
two combinations do not commute, so they cannot be considered as describing a single 
phase angle. Despite this deficiency, the Susskind-Glogower phase operators were widely 
used in description of optical fields until recently. The eigenkets of the Susskind-Glogower 
operator (2.5), given by 



generate the resolution of the identity, and, despite their nonorthogonality, they can be 
used to form the probability operator measure applied to the phase description by Shapiro 
and Shepard [1991]. 

Garrison and Wong [1970] made an attempt to construct a Hermitian phase oper- 
ator in the infinite- dimensional Hilbert space which, as they demanded, should satisfy 
the canonical commutation relation (2.1). Their work was almost completely forgotten. 
Bergou and Englert [1991] have made a comparison of the Garrison- Wong and Pegg- 
Barnett phase operators, pointing out that if the limit to infinite-dimensional space is 
performed on the latter operator (but not on the expectation values), the former operator 
is obtained. In their view the Pegg-Barnett phase formalism is only an approximation 
to the "correct" phase formalism. Bergou and Englert introduced their own quantum- 
phase description, which has not gained broader acceptance. Nevertheless, their paper is 
an essential contribution which clarifies a number of points in the quantum-phase prob- 
lem. The Garrison- Wong approach turned out to lead to phase distributions which are 
asymmetric and difficult to accept on physical grounds (Gantsog, Miranowicz and Tanas 
[1992]); e.g., even the vacuum is phase-anisotropic. For reference, we give a sketch of their 
approach in Appendix A. 

The renewed interest in quantum-phase problems has resulted in a reexamination of 
some of the earlier approaches and the creation of other, completely new descriptions. 
From a number of different phase formalisms that are now available, wc choose only two, 
which we shall apply for the description of the phase properties of fields generated in 
various nonlinear optical processes. These are: the Pegg-Barnett phase formahsm, which 
represents the canonical phase formalism based on the idea of finding a Hermitian phase 
operator, and another formalism based on the description of the optical phase through s- 
parametrized phase distributions, which for some values of s represents the experimentally 
measured phase probability distributions. 

2.1 Pegg-Barnett phase formalism 

Pegg and Barnett [1988, 1989] (and Barnett and Pegg [1989]) introduced the Hermitian 
phase formalism, which is based on the observation that, in a finite-dimensional state 
space, the states with well-defined phase exist (Loudon [1973]). Thus, they restrict the 
state space to a finite {a + 1) -dimensional Hilbert space H^^^ spanned by the number 
states |0), |1),..., |cr). In this space they define a complete orthonormal set of phase 



1 



oo 




(2.8) 
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states by: 



\0m) = —/== Y^QWi}nem)\n), m = 0, 1, . . . ,C7, (2.9) 
where the values of 6^ are given by 

Om = ^0 + ^. (2.10) 
cr + 1 

The value of 6*0 is arbitrary and defines a particular basis set of (cr+1) mutually orthogonal 
phase states. The number state \n) can be expanded in terms of the \6m) phase state basis 



as 



ri) = Y.\^rn){em\n) = ^==Y.^M-'^riem)\em). (2.11) 

71=0 V(7 + i 



m=0 



From eqs. (2.9) and (2.11) we see that a system in a number state is equally likely to be 
found in any state |^rn); and a system in a phase state is equally likely to be found in any 
number state \n). 

The Pegg-Barnett (PB) Hermitian phase operator is defined as 

= (*ioOpB = E^-l^-)(^-l- (2-12) 

m=0 

Of course, the phase states (2.9) are eigenstates of the phase operator (2.12) with the 
eigenvalues 9^, restricted to lie within a phase window between 6q and 6Q-\-2'Ka /{a-\-l). The 
Pegg-Barnett prescription is to evaluate any observable of interest in the finite basis (2.9), 
and only after that to take the limit cr — )> oo. 

Since the phase states (2.9) are orthonormal, {9m\0m') = ^mm', the fcth power of the 
Pegg-Barnett phase operator (2.12) can be written as 

= (2.13) 

m=0 

Substituting eqs. (2.9) and (2.10) into cq. (2.12) and performing summation over m yields 
explicitly the phase operator in the Fock basis: 

S> ^ I I 27r ^ exp[i(n-f?/)go] |n)(77/| 

' ° ^ a + 1 ^ a + l4;.exp[i(n-n0 27r/(a + l)]-l- > 

It is readily apparent that the Hermitian phase operator has well-defined matrix ele- 
ments in the number-state basis and docs not suffer from such problems as the original 
Dirac phase operator. A detailed analysis of the properties of the Hermitian phase oper- 
ator was given by Pcgg and Barnett |1989] and Barnett and Pcgg [1989]. 

The unitary phase operator exp(i$e) can be defined as the exponential function of the 
Hermitian operator ^q. This operator acting on the eigenstate 1^^) gives the eigenvalue 
exp(i^^), and can be written as (Pegg and Barnett [1988, 1989]): 

cr-l 

exp(il>e) = X l^)(^ + l|+exp[i((7 + l)^o]|(7)(0|, (2.15) 

n=0 
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Its Hermitian conjugate is 

[exp(il»e)]t = exp(-il'e) (2.16) 

with the same set of eigenstates \0m) but with eigenvalues exp(— i^^). This is the last term 
in eq. (2.15) that distinguishes the unitary phase operator from the Susskind-Glogower 
phase operator (Susskind and Glogower [1964]). The first sum in eq. (2.15) reproduces 
the Susskind-Glogower phase operator if the limit cr ^ oo is taken. In contrast to the 
Pegg-Barnett unitary phase operator, the Susskind-Glogower exponential operator (2.5) 
is defined as a whole and is not unitary, as is apparent from eq. (2.7). The sine and 
cosine operators in the Pegg-Barnett formalism are the sine and cosine functions of the 
Hermitian phase operator ^q. They are more consistent with the classical notion of phase 
than their counterparts in the Susskind-Glogower phase formalism. In particular, they 
satisfy the "natural" relations: 

cos^ $e + sin^ $e = 1, (2.17) 
[cosl>e,sin$e] = 0, (2.18) 
(n| cos^l>0|n) = (n| sin^l»0|n) = ^. (2-19) 

The last relation is also true for the vacuum state, in sharp contrast to the Susskind- 
Glogower phase operators. This is consistent with the phase of vacuum being random, as 
well as for any other number state. 

The Pegg-Barnett phase operator (2.14), expressed in the Fock basis, readily gives the 
phase-number commutator (Pegg and Barnett [1989]): 

ri 27r (n — n') exp[i(n — n')^ol , v , /, 

= Y — !l— ' L^\n)ln'\ 2.20 

^ ^ (7+l„^,exp[i(n-n0 27r/((7 + l)]-l' ' ^ > 

Equation (2.20) looks very different from the famous Dirac postulate of the phase-number 
commutator (2.1). Santhanam [1976] and Pegg, Vaccaro and Barnett [1990] examined 
canonically conjugate operators in the finite-dimensional Hilbert space. According to 
the generalized definition by Pegg, Vaccaro and Barnett, the photon number and phase 
are indeed canonically conjugate variables, similar to momentum and position or angular 
momentum and azimuthal phase angle. 

As the Hermitian phase operator is defined, one can calculate the expectation value 
and variance of this operator for a given state of the field |/). Moreover, the Pegg-Barnett 
phase formalism allows the introduction of the continuous-phase probability distribution 
which is a representation of the quantum state of the field and describes the phase prop- 
erties of the field in a very spectacular fashion. Examples of such phase distributions for 
particular states of the field will be given later. 

A general pure state of the field mode with a decomposition 

l/) = Ecn|n), (2.21) 

n=0 

can be re-expressed in the phase state basis, according to eq. (2.11), as 

I/) = ^:L=^^c„exp(-ine^)|ej. (2.22) 
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We should remark here that the coefficients c„ in the decomposition (2.21) in a finite- 
dimensional space should differ from the coefficients in the infinite-dimensional space if the 
state I/) is to be normalized. A short discussion of this problem is given in Appendix B. 
The phase probability distribution is given by (Pegg and Barnett [1989]): 



\{Om\f)f = 



1 



(7+1 



^c„exp(-in6'„ 



which leads to the expectation value and variance of the phase operator: 

m 
m 

If the field state |/) is a partial phase state, i.e., if the amplitude has the form 

the phase probabihty distribution (2.23) becomes 

1 



(2.23) 

(2.24) 
(2.25) 

(2.26) 



\{Omm - 



a+1 
1 



E6„exp[in((/? - 61^)] 



— — H — J2 ^nbk cos [{n -k){(p- 6'„ 



(2.27) 



The mean and variance of will depend on the chosen value of ^o- Judge [1964], in his 
description of the uncertainty relation for angle variables, suggested that the choice of 
phase window which minimizes the variance can be used to specify uniquely the mean 
or variance. For the partial phase state, the most convenient and physically transparent 
way of choosing is to symmetrize the phase window with respect to the phase (p. This 
means the choice 



and after introducing a new phase label 



a + 1 



the phase probability distribution (2.27) becomes 

1 2 

. 1 

n>k 



a 



a 



cos 



(n-k) 



271H 
a + 1 



(2.28) 



(2.29) 



(2.30) 



with II, which goes in integer steps from —a/2 to a/2. Since the distribution (2.30) is 
symmetrical in /i, we immediately get, according to eqs. (2.28)-(2.30): 



(/l<^>.|/) = ^• 



(2.31) 
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This result means that for a partial phase state with phase </?, the choice of ^ in 
eq. (2.28) relates directly the expectation value of the phase operator with the phase (p. 
With this choice of 60 the variance of the phase operator has a particularly simple form: 

ii^^ef) = E l^'mf)f. (2.32) 

^) ij,=-a/2 

So-called physical states play a significant role in the Pegg-Barnett formalism. Physical 
states \p) are defined by Pegg and Barnett [1989] as the states of finite energy. Most 
of the expressions in the Pegg-Barnett formalism take a much simpler form for physical 
states. For example, the commutator (2.20) reduces to 

[^e,n]p = -i[l-{a+l)\9o){9o\], (2.33) 

a form more similar to the standard canonical commutation relation (2.1). On the other 
hand, when physical states are considered, we can simplify the calculation of the sum in 
eqs. (2.24) and (2.25) by replacing it by the integral in the limit as a tends to infinity. 
Since the density of states is {a + l)/27r, we can write the expectation value of the kth 
power of the phase operator as 

9o+2n 

imelf) = / d9 9'P{9), (2.34) 
do 

where the continuous-phase distribution P{9) is introduced by 

P{9) ^ PMO) = jim (2.35) 

and 9m has been replaced by the continuous-phase variable 9. If the state |/) has the 
number-state decomposition (2.21), then the Pegg-Barnett phase distribution is (Pegg 
and Barnett [1989]): 



2Re E CmC^exp[— i(m — n)9] 



7i,m 

m'>n 



(2.36) 



In the case of fields being in mixed states described by the density matrix p, formula (2.36) 
generalizes to 

P{9) = ;i|l + 2Re E Pmnexp[-i(m-n)^]|, (2.37) 

where pmn = {m\p\n) are the density matrix elements in the number-state basis. Equa- 
tions (2.36) or (2.37) can be used for calculations of the Pegg-Barnett phase distribution 
for any state with known amplitudes c„ or matrix elements pmn- Despite the fact that 
they are exact, they can rarely be summed up into a closed form, and usually numerical 
summation must be performed to find the phase distribution. Such numerical summa- 
tions have been widely applied in studying the phase properties of optical fields. The 
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Pegg-Barnctt phase distribution, cqs. (2.36) or (2.37), is obviously 27r-pcriodic, and for 
all states with the density matrix diagonal in the number states, the phase distribution is 
uniform over the 27r-wide phase window. These are nondiagonal elements of the density 
matrix that lead to the structure of the phase distribution. The Pegg-Barnett distribution 
is positive definite and normalized. 

After introducing the continuous-phase distribution P{6), formula (2.32) for the phase 
variance, if the symmetrization is performed, can be rewritten into the form: 

TT 

{{A^ef) = J e^p{e)de. (2.38) 

— TT 

This means that as the phase distribution function P{0) is known, all quantum-mechanical 
phase characteristics can be calculated with this function in a classical-like manner. The 
result for the variance (2.38) is (Barnett and Pegg [1989]): 

TT^ (-1)'"--'' 

{{A^ef) = y + 4EU.|;^. (2.39) 

The value 7r^/3 is the variance for the uniformly distributed phase, as in the case of a 
single number state. 

For physical states there are some additional useful relations between the expectation 
values of the Pegg-Barnett phase operators and of the Susskind-Glogower phase operators. 
For example, the following relations hold (Vaccaro and Pegg [1989]): 





(cxp(im$sG))p, 


(2.40) 


{cos^e)p = 


(cOS$gG)p, 


(2.41) 


(sin l>e)p = 


(sin$sG)p, 


(2.42) 


(cos^$e)p = 


(c5s^$sg)p + ^((|0)(0|))„ 


(2.43) 


(sin2$e)p = 


(s'in'$sG)p + ^((|0)(0|))^. 


(2.44) 



where the subscript p refers to a physical state expectation value. 

2.2 Phase distributions from s-parametrized quasidistributions 

According to Cahill and Glauber [1969a, 1969b], the s-parametrized quasidistribution 
function W'^^^a) describing a field state, can be derived from the formula 

Vy(^)(a) = -Tr{pf^'\a)}, (2.45) 

TT 

where the operator T^^\a) is given by 

f'(^)(«) = 1 / exp« - a*0^^'HO d'e, (2.46) 

TT J 

L»(^)(^) = e^l«l'/2£)(^) (2.47) 
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with D{^) being the displacement operator; p is the density matrix of the field, and we 
have introduced the extra I/tt factor with respect to the original definition of Cahill and 
Glauber [1969b]. The operator T^^\a) can be rewritten in the form (Cahill and Glauber 
[1969a]): 

which gives explicitly its s-dependence. So, the s-parametrized quasidistribution function 
W*^*^ {a) has the following form in the number-state basis: 

W^^\a) = -^pmn{n\f^'\a)\m), (2.49) 

" m,n 

where the matrix elements of the operator (2.46) are given by (Cahill and Glauber [1969a]): 

X e-'^^^-'^^Vr^^'exp (^-p^^ ^n"" {t^) ' (2-^°^ 

in terms of the associate Laguerre polynomials LJ"^"(x). In eq. (2.50) we have also 
separated explicitly the phase of the complex number a by writing 

a = \a\e'^. (2.51) 

In the following, the phase 6 will be treated as the quantity representing the field phase. 

With the quasiprobabihty distributions W'^^^a) the expectation values of the s-ordered 
products of the creation and annihilation operators can be obtained by proper integrations 
in the complex a plane. In particular, for s = 1, 0, —1, the s-ordered products are normal, 
symmetrical, and antinormal ordered products of the creation and annihilation oper- 
ators, and the corresponding quasiprobabihty distributions are the Glauber-Sudarshan 
P-function, Wigner function, and Husimi Q-function. 

By virtue of the relation inverse to eq. (2.49), given by (Cahill and Glauber [1969b]): 

p = j d^af^-'\a)W^'\a), (2.52) 

the Pegg-Barnett phase distribution (2.37) can be related to the s-parametrized quasidis- 
tribution function (2.45) as follows (Eiselt and Risken [1991]): 

P(e) = j d^aK^'\a,e)W^'\a), (2.53) 

where the kernel is given by 

K^'\a,9) = — ^(n|r(-^)(a)|m)exp[i(m-n)^] (2.54) 

m,n 

in terms of the matrix elements (2.50) for (— s). The kernel (2.54) is convergent for 
s > — 1 only. Nevertheless, the remaining relation between the Husimi Q-function and 
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the Pegg-Barnett distribution can also be expressed by eq. (2.53), albeit with the following 
kernel (Miranowicz [1994]): 

n,m-' Vnlml 
X exp[i(n - m)9 - |q;|^ + a^* - q;*/3]. (2.55) 

On integrating the quasiprobability distribution W^^\a) over the "radial" variable \a\, 
we get the "phase distribution" associated with this quasiprobability distribution. The 
s-parametrized phase distribution is thus given by 

oo 

p('\e) = J W^'\a)\a\ d\a\ (2.56) 



or equivalently by 

oo 

pi') (9) = -jw^'\W,9)dW, (2.57) 





where integration is performed over the intensity W — \oi\^- Inserting eq. (2.49) into 
eq. (2.56) yields 

X / |ar""exp (^r^) l"l (2.58) 

If the definition of the Laguerre polynomial is invoked, the integrations in eq. (2.58) can 
be performed explicitly, and we get for the s-parametrized phase distribution a formula 
similar to the Pegg-Barnett phase distribution (2.37), which reads 



P'-'X^) = ;^|l + 2Re^p^,e-^("^-'^)^GW(m,n)|. 

L m>n J 



(2.59) 



The difference between eqs. (2.37) and (2.59) lies in the coefficients Gi^\m,n), which are 
given by 

/ 9 \ mm{m,n) , x i 



r(^-/ + i) 



The s-parametrized coefficients Gi'^\m,n), [eq. (2.60)], can be rewritten in a compact 
form (Miranowicz [1994], Leonhardt, Vaccaro, Bohmer and Paul [1995]) (m > n): 



n\ /s + 1\" / 2 



ml \s — ij \1 — s 
X r f ^ + l) '-(-+'^)) f ) (2.61) 



s+1 
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in terms of the Jacobi polynomials Pj^'t^) (x) , or equivalently as 
^(,), , 1 /s + iw 2 ^(— 




nl (m — n)\ \s — 1/ \1 — s 



X + ,F,[-n,^ + l,m-n+l,^) (2.62) 

using the hypergeometric (Gauss) function 2-Fi(a, c, x). 

For s=0, we have the coefficients for the Wigner phase distribution P^^\9), i.e., 
the phase distribution associated with the Wigner function. In this special case of 
s = 0, eq. (2.60) reduces to the expression obtained by Tanas, Murzakhmetov, Gantsog 
and Chizhov [1992], whereas eq. (2.62) goes over into the expression of Garraway and 
Knight [1992, 1993]: 

,,,, , f 2^^--y'M^YM^ rieven , , 

G^^'Hm.n) = { , ^/^Jr^l^L,, 2.63 

^ ' ^ I 2(m-n)/2 /^ r (rn+l / ) I ) 

V \ ml [(n— 1)/2]! 

Equations (2.61)-(2.63) are given for m > n. Otherwise the indices m and n should be 
interchanged. 

For s — —1, only the term with I — survives in eq. (2.60), and we get the coefficients 
for the Husimi phase distribution P^~^\9), i.e., the phase distribution associated with the 
Q-function. Now, eq. (2.60) reduces to (Paul [1974], Tanas, Gantsog, Miranowicz and 
Kiehch [1991], Tanas and Gantsog [1992b]): 

r f + 1) 

G^-^\m,n) = ^ g ^ (2.64) 

yn\m\ 

It is apparent from eqs. (2.59)-(2.62) that for the phase distribution associated with the 
P-function (>s=l) the coefficients G^^\m,n) become infinity, and one can conclude that 
the phase distribution P'^^^{9) is indeterminate. However, at least for a special class of 
states, summation can be performed numerically or even analytically for P^^\9). For 
instance, for the states described by the density matrix p of the form 

Pmn = |pmn|exp[-i(m - n)^?o], (2.65) 
the s-parametrized phase distribution P^^\9) can be rewritten as (Miranowicz [1994]): 

pW(^) = _ J 1 + 2 V cos[m(^ - t?o)] ^ (2.66) 

271- I m=l J 

with the coefficients 

oo 

'(^) = J2\Pm+n,n\G^'Km + n,n). (2.67) 



n=0 



Equations (2.66) and (2.67), for s = and ■j?o = 0, go over into expressions obtained 
by Bandilla and Ritze [1993]. Numerical calculation of lima^^ is usually straightforward. 
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Figure 1: The coefficients G*^'')(m, n) for (a) s = 0, and (b) s = —1 

For coherent states, the coefficients are equal to unity. Hence, PcoliO), given by 
eq. (2.66), is the Dirac delta function 6(9 — ■j?o) [see §3.1]. 

Formulas (2.59)~(2.62) allow calculation of the s-parametrized phase distributions 
for any state with known p^n said their comparison with the Pegg-Barnett phase dis- 
tribution, for which G^^\m,n) = 1. The phase distributions associated with particular 
quasiprobability distributions have been used widely in the literature to describe the phase 
properties of field states. For example, the Husimi phase distribution P'^~^\9) was used 
by Bandilla and Paul [1969], Paul [1974], Freyberger and Schleich [1993], Freyberger, Vo- 
gel and Schleich [1993a, 1993b], Leonhardt and Paul [1993a], Bandilla [1993] and KC94a 
in their schemes for phase measurement. Braunstein and Caves [1990] applied P^~^\9) to 
describe the phase properties of generalized squeezed states. The Wigner phase distribu- 
tion P^^\9) was used by Schleich, Horowicz and Varro [1989b, 1989a] in their description 
of the phase probability distribution for highly squeezed states. Herzog, Paul and Richter 
[1993] showed in general that the Wigner phase distribution can be interpreted as an 
approximation of the Pegg-Barnett distribution. To estimate the difference between the 
P'^'^\9) and P{9) they analyzed the deviation of the Wigner function W*^°(q;) for a phase 
state from Dirac's delta function. Recently, Hillery, Freyberger and Schleich [1995] have 
compared the Pegg-Barnett, Husimi, and Wigner phase distributions for large-amplitude 
classical states. Eiselt and Risken [1991] applied the s-parametrized quasiprobability dis- 
tributions to study properties of the Jayncs-Cummings model with cavity damping. 

For some field states the quasiprobability distribution functions >V*-*^(q;) can be found 
in a closed form via direct integrations according to the definitions (2.45)-(2.47), and 
sometimes the next integration leading to the s-parametrized phase distributions can 
also be performed according to definition (2.56). In the next Sections, we shall illus- 
trate the differences between the Pegg-Barnett phase distribution and the s-paramctrized 
phase distributions obtained by integrating the s-parametrized quasiprobability distribu- 
tion functions. For any field with known number-state matrix elements of the density 
matrix, the s-parametrized phase distribution can be calculated according to eq. (2.59) 
with the coefficients G''^^^(m,n) given by eq. (2.60). The distribution of the coefficients 
(m, n), for s = 0, —1, is illustrated in fig. 1. It is apparent that for the Husimi phase 
distribution the coefficients decrease monotonically as we go further away from the di- 
agonal. This means that all nondiagonal elements are weighted with numbers that 
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Figure 2: The minima of tlie Wigner phase distributions P^^\9), eq. (2.70), for the 
superpositions of two number states (2.68) for various values of n and m — n—1, 2, 3, 4. 



are less than unity, and the phase distribution for s = —1 is always broader than the 
Pegg-Barnett phase distribution (for which G^^\m,n) = 1). For s — the situation is 
not so simple, because the coefficients G^^\m,n) show even-odd oscillations with values 
that are both smaller and larger than unity. This usually leads to a phase structure 
sharper than the Pegg-Barnett distribution. Moreover, since the Wigner function (s = 0) 
can take negative values, the positive definiteness of the Wigner phase distribution is not 
guaranteed. Also, the oscillatory behavior of the coefficient G^°^(m, n) suggests that, at 
least for some states, the Wigner phase distribution P*^°^(^) can exhibit negative values. 
This nonclassical feature of P^'^\9) was shown explicitly by Garraway and Knight [1992, 
1993] for the simple example of the number-state superposition (only for convenience, we 
assume that m> n): 

I*) = 2-'^'^{\n) + \m)). (2.68) 

In a straightforward manner, the general expressions for the phase distributions P{0), 
[eq. (2.37)] and p(*)(^),[eq. (2.59)] reduce to 

P[d) = — (1 + cos[(m -n)^]), (2.69) 

and 

pW(^) = i- (l + G'(^)(m,n)cos[(m-n)^]) . (2.70) 

respectively. The Pegg-Barnett, P{0), and Husimi, P^~^\9), phase distributions are 
obviously positive definite for any superposition (2.68). As seen in fig. 2, the Wigner 
phase distribution P^'^\6) is positive for superpositions with odd m — n. However, it 
takes negative values for even m — n. In this case, the smaller is n for fixed m — n, or the 
higher is the value of m — n for a given n, the minimum of the Wigner phase distribution is 
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more strongly negative. Hence, one obtains the greatest negativity for the superposition 
(|0) + \2m))/y/2 in the hmit of m — )■ oo. As was emphasized by Garraway and Knight 
[1993] (see also fig. 2), for large values of n the Pegg-Barnett distribution is approached 
for both even and odd m. 

It is highly illustrative to consider analytically the special case of eq. (2.68) when 
m—n = 2 (Garraway and Knight [1992, 1993]). These rcsTilts will be helpful in the analysis 
presented in §3.2 for even and odd coherent states. Now, the coefficients G^^\m^ n), given 
by eqs. (2.60)-(2.62), can be rewritten in a much simpler form: 



1 - 



2^(n + l)(n + 2) 



s + 1 



n+2 



n + 1 



(2.71) 



For s = 0, eq. (2.71) goes over into (Garraway and Knight [1992, 1993]): 

'n + 2 
.n + 1 



G^^\n + 2,n) = ^) (2.72) 



and for s = — 1 one obtains 

G<-..(„ + 2,.) = (!i±i)"\ (2.73) 

Equation (2.72) provides direct proof of the oscillatory behavior of G^^\n + 2,n) with 
increasing n. For even n, the right-hand-side of eq. (2.72) is greater than unity, which 
implies a negative minimum of the Wigner phase distribution (2.70) [solid circles in fig. 2]. 
However, for odd n, the coefficients G^^\n+2, n) are less than unity and equal to G'^~^\n-\- 
2,n). Hence the Husimi and Wigner phase distributions for such states (with odd n) are 
equal and positive definite. 

From the form of the coefficients G'^^\m^n) it is evident that there is no s such that 
G^^\m,n) — 1 for all m, n. This means that there is no "phase ordering" of the field 
operators; that is, the ordering for which P^^\0) would be equal to P{0). However, for a 
given state of the field one can find s such that the two distributions are "almost identical" . 
Formula (2.59) is quite general, and it was used in earlier studies of the phase properties of 
the anharmonic oscillator (Tanas, Gantsog, Miranowicz and Kielich [1991]), parametric 
down conversion (Tanas and Gantsog [1992b])), and displaced-number states (Tanas, 
Murzakhmetov, Gantsog and Chizhov [1992]). A disadvantage of formula (2.59) is the 
fact that the numerical summations can be time consuming and even difficult to perform 
for field states with slowly converging number-state expansions. This, for example, is 
the case for highly squeezed states. In some cases, instead of using the number-state 
expansions, we can find analytical formulas for P^^\d) via direct integrations, as shown 
in §3. In many cases such formulas can be treated as good approximations to the Pegg- 
Barnett phase distribution. 



3 Phase properties of single-mode optical fields 

Optical fields produced in nonlinear optical processes have specific phase properties which 
depend on the nonlinear process in which the field is produced and on the state of the field 
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before it undergoes the nonlinear transformation. Since there is a large variety of nonlinear 
optical processes, there is the possibility to generate fields with different phase properties. 
Here, some examples of such field states will be given and their phase properties discussed 
briefly. 



3.1 Coherent states 

The most common single-mode field in quantum optics is a Glauber coherent state. Its 
phase properties have probably been analyzed within each known phase formalism. We 
shall focus our attention on two of them only. 

The s-parametrized quasiprobability distribution function for a coherent state 

|ao) = ^(«o)|0) (3.1) 

can be calculated from eqs. (2.45)-(2.47) as 

= l|exp[(a-ao)r-(a*-aS)e + «ier/2](0|I)(e)|0)d2e 
■ J cxp[(a - ao)C - («* - ao)C + - |ao|V2] d'i 



1 



7r2 



1 2 r 2 , ,2i , , 

= exp<^-- \a-ao\ >. (3.2) 

The corresponding s-parametrized phase distribution is (Tanas, Miranowicz and Gantsog 
[1993]; for s = see also Garraway and Knight [1993] and BR93): 

oo 

p('\e) - J W^'\a)\a\ d\a\ 



^ exp[-(Xo^ - X')] {exp(-X2) + V^X{1 + erf (X))} , (3.3) 



where 



X = X^'\e) = y'-A-|ao|cos(^-i?o), 



(3.4) 



and Xq — X(*)('j?o); '&o is the phase of ckq. The phase distribution P^^\9) associated with 
the P-function can be obtained from eqs. (3.3) and (3.4) in the limit of s — > 1: 

p(i)(^) ^ 5(^-^o), (3.5) 

which is the Dirac delta function. This result can also be achieved from eq. (2.66). As 
was shown by Miranowicz [1994], the coefficients are unity for arbitrary m. Hence, 
eq. (2.66) reduces to 

1 f °° 1 

P^'H^) = ^ 1 1 + 2 E cos[m(^ - ^o)] [ 

= 5(^-^o), (3.6) 
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Figure 3: Phase distributions for the coherent states with the mean number of photons: 
(a) |q;oP = 2, and (b) |q;oP = 0.01; the Pegg-Barnett distribution (sohd hue), the Wigner 
phase function P^^\9) (dashed hue), and the Husimi phase distribution P^~^\9) (dotted- 
dashed hue). 



which is the desired function (3.5). This example shows that the general expression (2.59) 
for the s-parametrized phase distributions is also valid in the special case of s=l. 

Formula (3.3) is exact, it is 27r-periodic, positive definite and normalized, so it satisfies 
all requirements for the phase distribution. Moreover, formula (3.3) has a quite simple 
and transparent structure. For small |q;o|, the first term in braces plays an essential role, 
and for |q;o| ^ we get a uniform phase distribution. For large |ao|; the second term in 
the braces predominates, and if we replace erf(X) by unity, we obtain the approximate 
asymptotic formula given by Schleich, Dowling, Horowicz and Varro [1990] (for s = 0): 

P(°)(^) ^ ^|ao|cos(^-i?o)exp[-2|ao|^sin2(^-i?o)], (3.7) 

which however, can be applied only to — 7r/2 < (^ — ■J?o) < 7r/2. After linearization of 
formula (3.7) with respect to 9, the approximate formula for coherent states with large 
mean number of photons obtained by Barnett and Pcgg [1989] is recovered. The presence 
of the error function in cq. (3.3) handles properly the phase behavior in the whole range 
of phase values —n < {9 — i^o) < 

The Pegg-Barnett distribution P{9) for the coherent state \ao) can be calculated from 
eq. (2.36) with the superposition coefficients 

Cn = 6„exp(in^?o), K = exp(-|Q;or/2)— (3.8) 

vm 

The exact formula for the s-parametrized phase distributions P^'^\9) for coherent states 
is given by eqs. (3.3) and (3.4). Alternatively, the P^^\9) are given by eq. (2.59) after 
insertion of c„ given by eq. (3.8). In fig. 3 we show the phase distributions P{9), P'^'^\9), 
and P^~^\9) for a coherent state with the mean number of photons |q;oP = 2 (a), and 
I Clop = 0.01 (b). It is seen that the Pegg-Barnett phase distribution is located somewhere 
between the Wigner and Husimi phase distributions. It becomes closer to P(°)(^) for 
Ictop ^ 1, and closer to P^~^\9) for |aoP <^ 1. For |q;oP — )■ oo, the Pegg-Barnett dis- 
tribution tends to the Wigner phase distribution (Schleich, Horowicz and Varro [1989a], 
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Barnctt and Pcgg [1989]), and for |aoP all the distributions tend to the uniform 
distribution, but the Pegg-Barnett distribution in this region tends to the Husimi phase 
distribution. This means that for coherent states with large mean numbers of photons, 
P^^\9) is a good approximation to the Pegg-Barnett phase distribution, while for small 
numbers of photons P^~^\9) becomes a good approximation to the Pegg-Barnett distri- 
bution. 

3.2 Superpositions of coherent states 

Superpositions of macroscopically distinguishable coherent states have attracted much 
interest (see, for example, Buzek and Knight [1995] and references therein) due to their 
property of being prototypes for the Schrodinger cats, and important nonclassical proper- 
ties, such as sub-Poissonian photon statistics, quadrature squeezing, higher-order squeez- 
ing, etc.. Their phase properties have also been a subject of interest. 

Let us consider the normalized superposition l'^) of coherent states defined as 



This superposition of two well-separated components is called a Schrodinger cat, whereas 
for N > 2 the notions Schrodinger cat-like state or kitten states are often used. The phase 
distributions P(^), [eq. (2.37)] and P^'\e), [eq. (2.59)] for the state (3.9) can be rewritten 
in a form showing explicitly the superposition structure (Tanas, Gantsog, Miranowicz 
and Kiclich [1991], Garraway and Knight [1992, 1993], Buzek, Gantsog and Kim [1993], 
Buzek, Kim and Gantsog [1993], Tara, Agarwal and Chaturvedi [1993], Hach III and 
Gerry [1993], Buzek [1993], Miranowicz [1994], Buzek and Knight [1995]). 

The Pegg-Barnett phase distribution P{9) sphts into two terms (Tanas, Gantsog, 
Miranowicz and Kielich [1991]): 



N 




(3.9) 



p{e) 



Po{0) + Pint(^), 



(3.10) 



where 



JV 



^o(^) 



Elcfcl'n(^) 




k=l 



is the sum of phase distributions 




(3.12) 



for the coherent states of the superposition, and the second distribution 



N 



Pint{0) 



E CkC*iPki{e), 



(3.13) 



k,l=l 
k^l 



Pki{0) = ^Y.bmbne^p[im{(Pk + ^o-0)-in{(l>i + i9o-0)] 



(3.14) 



m,n 
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represents interference terms emerging from the quantum interference between the com- 
ponent states of the superposition. In fig. 4, the phase distributions (3.10) and (3.11) 
are presented in polar coordinates for the discrete superpositions of coherent states in 
the anharmonic oscillator model [see §3.5,eq. (3.57)]. It is evident from fig. 4 that as 
the number of components in the superposition becomes larger, the interference terms 
play an increasing role and the symmetry of the Pegg-Barnett distribution [eq. (3.10)] is 
destroyed. These terms are negligible for well-separated components of the superposition 
only (Tanas, Gantsog, Miranowicz and Kielich [1991]). 

Analogously, the s-parametrized quasidistribution W^^\a) for the superposition 
state (3.9) is represented as (Miranowicz [1994]; for s = — 1 see Miranowicz, Tanas and 
Kiehch [1990], and for s = see also Buzek and Knight [1995]): 



where the sum of coherent terms is 



N 



k=l 



with 



1 2 



TTl 



The interference part is given by: 



>VS(«) = 2 Y: \ck\ \ci\ Wl'i\a) 



exp 



N 



is). 



k,l=l 
k>l 



with 



1 2 

TT 1 — S 



exp 



1-s 



a 



+ 



a. - c'^'ao 



+ 2 



X cos|7fe 

4 , 



7z 



1 + g, 
1-s 

1 + s 



\ao\ sm 



1 



+ 



1-s 



\a\ \ao\ cos 



Q;o|^sin(0fe - 0^) 

k + 4>i 



+ -00- 9] sin 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



In eq. (3.19) the phases are 7^ = Argc^, 9 = Arga, 'do = Argao, and appears in the 
definition (3.9). On integration, we obtain the following form of the s-parametrized phase 
distribution P^'\9) (Miranowicz [1994]): 



P^^\9) = P^'\9) + Pt;{9), 



i.e., a simple sum 



N 



(3.20) 



(3.21) 



k=l 
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Figure 4: The Pegg-Barnett phase distributions in polar coordinates for the discrete 
superpositions of coherent states (3.57) with N—2-7 components in the anharmonic os- 
cillator model for the initial mean photon number |q;oP = 4; the exact phase distributions 
P{9), cq. (3.10) [sohd lines], and the functions Po{0), eq. (3.11) [dashed lines], without 
interference contribution. 

of coherent terms 

Pk^'\0) = ^ exp[-(X2 - Xl)] {exp(-X|) + V^X^ [1 + erf(X,)]} , (3.22) 

where 

cos(^^-^o-0fc), Xo = X(j(i?o + 0fe) (3.23) 

= Ec^^rn^fV) (3.24) 

{X^Kkl - Xli)] 

-Xl) + ^Xki [1 + erf(XfcO]} , (3.25) 
with 

X^i^X^^iO) = iyiII|ao|{exp[i(0, + ^o-^)] 

+ exp[-i((/), + ^?o-^)]}, (3.26) 



and the sum 

of the interference terms 

= ^exp[- 
X |exp(- 
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1 — s 1 + s 

«w = + exp{i(0fe - 0;)}, (3.27) 

The Schrodinger cat of the form: 

IV^) = |a,7) = AA^(|a)+exp(i7)|-a)), (3.28) 
with the normahzation 

A/; = [2 (l + cos7exp(-2|aH)] , (3.29) 

is a special case of the superposition state (3.9). The cat (3.28) consists of two coherent 
states la) and \ — a), which are = tt out of phase with respect to each other. The 
state (3.28) is not only of theoretical interest, since several methods were proposed for 
generation and measurement of this Schrodinger cat (e.g., Brune, Haroche, Raimond, 
Davidovich and Zagury [1992], Garraway, Sherman, Moya-Cessa, Knight and Kurizki 
[1994]). The state \a,7i/2) (i.e., for 7 = 7r/2) is called the Yurke-Stoler coherent state 
(Yurke and Stoler [1986, 1988]). This state can be generated in the anharmonic oscillator 
model (see §3.5). For other choices of 7, the state (3.28) goes over into the even coherent 
state \a,0) or the odd coherent state \a,n), which have the following Fock representa- 
tions (Pefina [1991]): 

00 2n 

\a,0) = cosh-^/'(|any ^=|2n), (3.30) 

n=o J{2n)\ 



00 ^2n+l 

\a,7T) = sinh-V'(|a|^) V " |2n + l). (3.31) 

n=0 ^^(271+1)! 

The dissimilar phase properties of the even/odd coherent states were analyzed by Gar- 
raway and Knight [1993] (see also Buzek and Knight [1995]). Their phase distributions 
P{9) and P^'\9) can be obtained readily from the general expressions (3.10) and (3.20), 
respectively. Obviously, they consist of the normalized sum of the phase distributions 
Pi,2{d) (or Pi^li^^)) coherent states located at a and {—a) in the phase space plane, 
together with an additional interference term Puid) (or Pi2{(^))- As was shown by Gar- 
raway and Knight [1993], the Wigner phase distribution P(°)(6') for the even coherent 
state [eq. (3.30)] can exhibit negative values, in contrast to P^^\9) for the odd coherent 
state [eq. (3.31)], which never does. The Fock expansion [eq. (3.30)] of the even coherent 
state contains only even photon numbers similar to the superposition |ncvcn) + l^cvcn + 2) 
discussed by us in §2.1 [see eq. (2.68) and fig. 2]. Analogously, the odd coherent state 
[eq. (3.31)] and |nodd) + l^odd + 2) contain only odd number states. Hence, these dissimi- 
lar features of the functions P(°)(^) for \a, 0) and \a, n) are well understood for the same 
reasons as those given in §2.1 in the analysis of the Wigner phase distribution for the 
superposition of the two number states and the interpretation of the oscillatory behavior 
of the coefficients G^^\m,n) [fig. la]. 
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3.3 Squeezed states 

Squeezed states have phase-sensitive noise properties, and it is particularly interesting to 
study their phase properties. Sanders, Barnctt and Knight [1986], Yao [1987], Loudon and 
Knight [1987], and Fan and Zaidi [1988] used the Susskind-Glogower formalism in a de- 
scription of the phase fluctuations of squeezed states. Lynch [1987] applied the measured- 
phase formalism of Barnctt and Pcgg [1986]. Vaccaro and Pegg [1989] and Vaccaro, 
Barnett and Pegg [1992] investigated phase properties of a single-mode squeezed state 
from the point of view of the new Pegg-Barnett phase formalism. Gr0nbech- Jensen, Chris- 
tiansen and Ramanujam [1989] made comparisons of the phase properties of a single-mode 
squeezed state obtained according to different phase formalisms, including that of Pcgg 
and Barnett. Burak and Wodkiewicz [1992] introduced a phase-space propensity descrip- 
tion of quantum-phase fluctuations and analyzed, in particular, squeezed vacuum. The 
phase properties of the squeezed states have recently been studied by Cohen, Ben-Aryeh 
and Mann [1992], and by Collett [1993a, 1993b]. Various measures of phase uncertainty 
and their dependence on the average number of photons were studied by Freyberger and 
Schleich [1994]. 

Squeezed states (ideal squeezed states, two-photon coherent states) are defined by 
(see Loudon and Knight [1987]): 

|«o,C) = ^(ao)^(C)|0), (3.32) 
where S{Q is the squeezing operator 

•^(C) = exp Qca^ - ICS^^) , (3.33) 
and ( is the complex squeeze parameter 

C = re''\ r = |C|. (3.34) 

The direct integrations lead to the s-parametrized quasiprobability distribution (for rj — 
0): 



f 2 2 1 

X exp < [Im(Q; — cto)]^ [Re(a — cto)]^ \ , (3.35) 

[ H-s 1 -s J 

where we have used the notation /i = e^^ . After integration over assuming that cto is 
real, we arrive at the formula (Tanas, Miranowicz and Gantsog [1993]): 



^ ^ 27r (/i - s) cos2 9 + (/x-i - s) sin^ 9 

X exp[-(X2 - X^)] {exp(-X2) + ^X (1 + erf (X))} , (3.36) 

where 



X = Xi^\9) = J^- -oVJ^scos9 

V /^"' - s ^{i_i-s) cos2^ + (At-i-s) sin^e 
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Although the variable X is slightly different, the main structure of the phase distribution 
is preserved. Formula (3.36) is valid for both small and large ctQ. For oq = we have 
the result for squeezed vacuum. After appropriate approximations, one easily obtains the 
formula derived by Schleich, Horowicz and Varro [1989a] for a highly squeezed state. 

The exact analytical formula for the s-parametrized phase distribution for squeezed 
states, as given by eqs. (3.36) and (3.37), for the squeezed vacuum takes the form: 

1 ^(// - s)(//-i - s) 

^'^^^ " 2^(//-s)cos2^ + (//-i-s)sin2^' ^^-^^^ 

where /i — exp(2r). This formula exhibits a two-peak structure with peaks for 9 — ±7r/2 
(for r > 0) . It is easy to find that the peak heights are 



meaning that for s — 0, the peak height is proportional to /i. One can easily check that 
the Pegg-Barnett result lies between the s = and s = — 1 results. Qualitatively, all three 
distributions give the same two-peak phase distributions, but they differ quantitatively: 
the sharpest peaks are those of P(°)(6'), and the broadest those of P^-^\9). 

For squeezed states with non-zero displacement cto, an additional factor of a form 
identical with that for coherent states, except for the different meaning of X{9), appears 
in the phase distribution P^^\9). Since this extra factor shows a peak at ^ = 0, a 
competition arises between the two-peak structure of the squeezed vacuum and the one 
peak structure of the coherent component. This competition leads to the bifurcation in 
the phase distribution discussed by Schleich, Horowicz and Varro [1989b, 1989a]. Figure 5 
illustrates such a bifurcation for ao — 1, bs exhibited by the Wigner and Husimi phase 
distributions plotted on the same scale to visualize the differences. Qualitatively, the 
pictures are quite similar, and differ only in the widths of the peaks. The Pegg-Barnett 
distribution in this case is very close to the Wigner phase distribution, and for this reason 
we omit it here. To calculate the Pegg-Barnett phase distribution one can apply formula 
(2.36) with c„ given by (see Loudon and Knight [1987]): 

ao + cto®^"' tanh r 



{n\aoX) 



\Jn\ coshr 





n/2 


^e^''^ tanh r 


Hn 


.2 





V2e2''? tanh r 

X exp|-^[|aor + aoV'''tanhr]|, (3.40) 

assuming 1] = (results for = 7r/2 can be obtained on replacing r by — r). 

Approximate analytical formulas for the phase variance as well as cosine and sine 
variances were obtained by Vaccaro and Pegg [1989] for weakly squeezed vacuum. For 
large squeezing the squeezed vacuum phase variance asymptotically approaches 7r^/4, 
which corresponds to the phase distribution with two symmetrically-placed delta- functions 

P{9) = \[S{9-7:/2) + 6{9 + 7:/2)]. (3.41) 

Ideal squeezed vacuum is generated in the parametric down-conversion process, in which 
the pump field is treated as a constant classical quantity. Taking into account the quantum 
character of the pump one finds that the signal field is no longer the ideal squeezed vacuum 
and its phase properties are different (Tanas and Gantsog [1992a]) [see §4.5]. 
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Figure 5: Pictures of the phase bifurcation for the squeezed state with the mean number of 
photons |q;oP = 1- The distributions are: (a) P^^\9), and (b) P^~^\9). The Pegg-Barnett 
distribution is very close to (a). 




Figure 6: The Pegg-Barnett phase distribution (3.44) of the Jaynes-Cummings model as 
a function of scaled time T = gt / {27:\ao\) for the initial mean photon number |q;oP = 20. 

3.4 Jaynes-Cummings model 

The Jaynes-Cummings model (Jaynes and Cummings [1963]) (see reviews by Yoo and 
Eberly [1985] and Shore and Knight [1993]) is the most popular model used to describe 
the resonant interaction of single two-level atom with single mode of the electromagnetic 
field. One of the most remarkable effects predicted theoretically (Eberly, Narozhny and 
Sanchez- Mondragon [1980], Narozhny, Sanchez- Mondragon and Eberly [1981]) and then 
observed experimentally (Rempe, Walther and Klein [1987]) in the Jaynes-Cummings 
model are collapses and revivals of the atomic inversion. Eiselt and Risken [1989], using 
the Q-function, have shown that the collapses and revivals can be understood in terms 
of interferences in phase space. Phoenix and Knight [1990] mentioned the splitting of the 
phase probability distribution into two counter-rotating satellite distributions in a model 
consisting of two degenerate atomic levels, coupled through a virtual level by a Raman- 
type transition. Dung, Tanas and Shumovsky [1990] discussed the collapses and revivals 
in this model from the point of view of the field-mode phase properties studied in the 
framework of the Pegg-Barnett formalism. 
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The model is described by the Hamiltonian (at exact resonance): 

H = huj{a^a + R^) + hg{R^a + Ra''), (3.42) 

where a"'" and a are the creation and annihilation operators for the field mode; the two- 
level atom is described by the raising, R\ and lowering, R, operators and the inversion 
operator R^, and g is the coupling constant. 

To study the phase properties of the field mode, we must know the state evolution of 
the system. After dropping the free evolution terms, which change the phase in a trivial 
way, and assuming that the atom is initially in its ground state and that the field is in a 
coherent state \ao), the state of the system is found to be 

oo 

li^it)) = ^Kewi^n-do) cos{^/ngt)\n,g) -ism{^/ngt)\n- l,e) , (3.43) 

n=0 

where \g) and |e) denote the ground and excited states of the atom, the coefficients 6„ are 
given by eq. (3.8) and -i^o is the coherent state phase. 

According to the Pegg-Barnett formalism, one obtains the phase distribution P{9) in 
the form (Dung, Tanas and Shumovsky [1990]): 

P{e) = — < 1 + 2 V bnbk cos[(n - k)9] cos[(v^ - Vk)gt] 1 . (3.44) 

[ n>k J 

This formula can be rewritten into the form: 

m = l[p+{o)+p-m, (3.45) 

where 

P±{e) = ^ll + 2Y,bnbkC0s[in-k)e^iV^-Vk)gt]\, (3.46) 

n>k ) 

which shows explicitly that as time elapses, the phase distribution P{9) splits into two 
distinct, right and left rotating, distributions in the polar coordinate system. Polar plots 
of the phase distribution are shown in fig. 6 (the time T = gt/{2-K\ao\) is scaled in the 
revival times). So, after a certain interval of time, the two counter-rotating distributions 
"collide", and at that time the components of the field oscillate in phase and one can 
expect the revival of the atomic inversion. The numerical calculations corroborate this 
statement (Dung, Tanas and Shumovsky [1990]). This behavior of the phase distribution 
resembles the behavior of the Q-function studied by Eiselt and Risken [1991]. The time 
behavior of the phase variance together with the phase probability density distribution 
carries certain information about the collapses and revivals. To show this, we first give 
the exphcit expression for the variance. Using eqs. (2.36) and (3.44), one obtains 

7]-2 ( — -Wn-k 

((A$,)2) = +4^^^UfeCos[(V^-V^)^t]. (3.47) 

n>k ~ 

Variance (3.47) is illustrated graphically in fig. 7 for |q;oP = 20. The variance goes up 
initially and reaches a maximum at the scaled time T — 1, which is the first revival time. 
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Figure 7: Evolution of (a) the mean photon number (n) and (b) the variance ((A$5i)^) of 
the Pegg-Barnett phase operator for the Jaynes-Cummings model as a function of scaled 
time T = gt/{2n\ao\) for |q;oP = 20. 

The revival times correspond to the extrema of the phase variance. In this way, the well 
known phenomenon of collapses and revivals has obtained clear interpretation in terms 
of the cavity-mode phase. More details can be found in the paper by Dung, Tanas and 
Shumovsky [1990]. The dynamical properties of the field phase in the Jaynes-Cummings 
model were studied by Dung, Tanas and Shumovsky [1991a], and the effects of cavity 
damping by Dung and Shumovsky [1992]. Some generalizations of this simple model 
were also considered from the point of view of their phase properties (Dung, Tanas and 
Shumovsky [1991b], Mcng and Chai [1991], Meystre, Slosser and Wilkens [1991], Dung, 
Huyen and Shumovsky [1992], Meng, Chai and Zhang [1992], Peng and Li [1992], Peng, 
Li and Zhou [1992], Wagner, Brecha, Schenzle and Walther [1992, 1993], Fan [1993], Jex, 
Matsuoka and Koashi [1993], Drobny, Gantsog and Jex [1994], Fan and Wang [1994], 
Meng, Guo and Xing [1994]). 

3.5 Anharmonic oscillator model 

The anharmonic oscillator model is described by the Hamiltonian 

H = huja)a + ^hKa^^a'^, (3.48) 

where a and are the annihilation and creation operators of the field mode, and k is 
the coupling constant, which is real and can be related to the nonlinear susceptibility x*^^^ 
of the medium if the anharmonic oscillator is used to describe propagation of laser light 
(with right or left circular polarization) in a nonlinear Kerr medium. If the state of the 
incoming beam is a coherent state |q;o), the resulting state of the outgoing beam is given 
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by 



IV'(t)) = U{T)\ao) 



= exp(-|ao|V2) ^exp 
n=o V n! 



i—nin 
. 2 ^ 



1) 



\n) 



(3.49) 



where r = —nt. In the problem of hght propagating in a Kerr medium, one can make the 
replacement t — —z/v to describe the spatial evolution of the field instead of the time 
evolution. On introducing the notation cto = |q;o| exp(i'j?o) the state (3.49) can be written 
as 



= l^^nexpii 



n=0 



w&o + —n{n — 1) 



\n) 



(3.50) 



where bn is given by eq. (3.8). 

The appearance of the nonlinear phase factor in the state (3.50) modifies essentially 
the properties of the field represented by such a state with respect to the initial coherent 
state |ao)- It was shown by Tanas [1984] that a high degree of squeezing can be obtained 
in the anharmonic oscillator model. Squeezing in the same process was later considered 
by Kitagawa and Yamamoto [1986], who used the name crescent squeezing because of the 
crescent shape of the quasiprobability distribution contours obtained in the process. The 
evolution of the quasiprobability distribution Q{a, a*) in the anharmonic oscillator model 
was considered by Milburn [1986], Milburn and Holmes [1986], Perinova and Luks [1988, 
1990], Daniel and Milburn [1989], and Perinova, Luks and Karska [1990]. The states 
that differ from coherent states by extra phase factors, as in eq. (3.50), are the general- 
ized coherent states introduced by Titulaer and Glauber [1966] and discussed by Stoler 
[1971]. Bialynicka-Birula [1968] has shown that, under appropriate periodic conditions, 
such states become discrete superpositions of coherent states. Yurke and Stoler [1986], and 
Tombesi and Mecozzi [1987] discussed the possibility of generating quantum-mechanical 
superpositions of macroscopically distinguishable states in the course of evolution for the 
anharmonic oscillator. Miranowicz, Tanas and Kiehch [1990] have shown that superposi- 
tions with not only even but also odd numbers of components can be obtained. 

The Pegg-Barnett Hermitian phase formalism has been applied to the study of the 
phase properties of the states (3.50) by Gerry [1990], who discussed the limiting cases 
of very low and very high light intensities, and by Gantsog and Tanas [1991f], who gave 
a more systematic discussion of the exact results. Phase fiuctuations in the anharmonic 
oscillator model were also analyzed within former phase formalisms (Gerry [1987], Lynch 
[1987]). 

The continuous Pegg-Barnett phase probability distribution (2.36) for the field 
state (3.50) takes the following form: 



PiO) = —\l + 2j2bnbkCos 

I n>k 



T 



(n - k)e - -[n{n - 1) - k{k - 1)] 



(3.51) 



and the s-parametrized quasiprobability distribution function is now given by (Miranow- 
icz [1994]): 



1-s 



exp 



1-s 



{\ao\'^ + \af) 
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X 



E 

m>n 



|a|™-"|aor+™ 



m! 



1-s 



s + 1 
s-1 



A\a\ 



X COS |(m — n)('j?o — ^) + ^[ m(m — 1) — n{n — 1) ]| 

4 \ 

+ Jo(iY3^|Qi| kol) J , 



(3.52) 



where Jo(2;) is the Bessel function. For r = 0, >V*^''-'(q;), given by eq. (3.52), is the coherent- 
state distribution [eq. (3.2)]. In the special case, for Q-function (s = —1), eq. (3.52) 
reduces to 

2 



Q{a,T) 



exp(— |Q!| 



E 

n=0 



exp ^i— n(n — 1) 



(3.53) 



The s-parametrizcd phase distribution, resulting from eq. (3.52) is 

L m>n 



X COS 



(m — n)(^ — -do) (m(m — 1) — ^1(71 — 1)) 



(3.54) 



where the coefficients G^^\m,n) are given by eq. (2.60). Symmetrization of the phase 
window with respect to the phase -do as done for the Pegg-Barnett phase distribution 
[eq.(3.51)] is equivalent to introduction of the relative phase variable 9 — -§0, and the two 
formulas differ only by the coefficients G^^\m,n), as in cq. (2.59). For r = 0, eqs. (3.51) 
and (3.54) describe the phase probability distributions for the initial coherent state |ao)- 
When the nonlinear evolution is on (r 7^ 0), the distributions P{9) and P^'\9) acquire 
some new and very interesting features. A systematic discussion of the properties as well 
as the plots of P{9) and P^'^\9) are given by (Tanas, Gantsog, Miranowicz and Kielich 
[1991]) and by (Gantsog and Tanas [19911]). 

The phase distribution P{9) can be used to calculate the mean and the variance of the 
phase operator, defined by eqs. (2.24) and (2.25). The results are (Gantsog and Tanas 
[1991^): 



TT 

{^|J(r)\^0\^|J(r))^'&o + J 9P(9)d9 



\n—k 



- 2 ^ hnh ^ ^^'\\ in\Un{n - 1) - k{k - 1)] \, 
n>k n- K 



(3.55) 



((Ad,)^ 



9^P{9) d9 - 



9P{9) d9 



Ir-TT 

\n—k 



n>k 



2EMfe- 

n>k 



^ sin |-[n(n — 1) — A;(A; — 1) 



(3.56) 
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Figure 8: Evolution of (a) the mean phase (3.55) and (b) the phase variance (3.56) for 
the anharmonic oscillator model. 



For r = 0, we recover the results for a coherent state with the phase 'do [eqs. (2.31) 
and (2.39)]. The nonlinear evolution of the system leads to a nonlinear shift of the mean 
phase and essentially modifies the variance. An example is illustrated in fig. 8, where 
the evolution of the mean phase (fig. 8a) and its variance (fig. 8b) are plotted against 
T for various values of Iccop- We have assumed 'j?o = 0, and the window of the phase 
values is taken between — vr and vr. The evolution is periodic with the period 27r, so the 
initial values are restored for r = 27r. Figure 8a shows the intensity-dependent phase 
shift. The amplitude of the mean phase oscillation becomes larger with increasing mean 
number of photons. The hne 7r^/3 in fig. 8b marks the variance for the state with random 
distribution of phase. It is seen clearly from fig. 8b that the stronger the initial field, the 
higher the phase variance. For |q;oP = 4 the phase variance increases rapidly and most of 
the period oscillates around 7r^/3 - the value for uniform phase distribution. This means 
that the phase is randomized during the evolution, although it periodically reproduces its 
initial values. This tendency is even more pronounced when the mean number of photons 
increases. The periodicity of the evolution is destroyed by damping (Gantsog and Tanas 
[1991b]). The sine and cosine functions of the phase were also calculated and the results 
compared with other approaches (Gantsog and Tanas [19911]). 

The local minima in ((A^g)^) apparent in fig. 8 indicate the points in the evolution 
in which superpositions of coherent states are formed, and the phase variance decreases 
at these points. This occurs for r = 27rM/N {N—2, 3, 4,. . . , and M, are mutually 
prime numbers), for which the P{0) and P^^\9) plotted in polar coordinates show A^-fold 
symmetry, confirming the generation of discrete superpositions of coherent states with 2, 
3, 4, . . . components: 



^(r = 2.-) 



27V 



where the phases (pk are simply 



TT 



J^Ckl exp(i0fc)tto), 

k=l 



A; = 1, . . . , 2A^, 



(3.57) 



(3.58) 
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and the superposition coefficients c^, representing the so-called fractional revivals, are 
given by (Averbukh and Perelman [1989], Tanas, Gantsog, Miranowicz and Kielich [1991]): 

1 / TT \ 

Cfe = ^12^ew[-^j;^[^k-Mn{n-l)]j. (3.59) 

Such superpositions, created during the evolution of the anharmonic oscillator model, 
have very specific phase properties discussed in §3.2. Plots the phase distributions (3.10) 
and (3.11) (where N should be replaced by 2A'") for the superpositions (3.58) with several 
components are presented in fig. 4. The phase distribution indicates the superpositions 
in a very spectacular way, as shown by Tanas, Gantsog, Miranowicz and Kielich [1991], 
Gantsog and Tanas [1991f] and Sanders [1992] for the anharmonic oscillator model, and 
by Paprzycka and Tanas [1992] for the model with higher nonlinearities. 



3.6 Displaced number states 

Other states which are interesting from the point of view of their phase properties are 
the displaced number states \ao,nQ) generated by the action of the displacement operator 
D{ao) on a Fock state \no) (see de Ohveira, Kim, Knight and Buzek [1990]); 

\ao,no) = D{ao)\no). (3.60) 

In a special case, when no — 0, the states (3.60) become a coherent state \ao). The 
s-parametrized quasiprobability distribution for the state (3.60) is 

>V(^)(a) = i^(_i)^>^l + ^^'^ 



ttI — s Vl — s 



whereas the s-parametrized phase distribution becomes (Tanas, Miranowicz and Gantsog 
[1993]): 



(ey \ n n 



k=o k\ \ 2 J \k 
xi:fn%^«-^^m-.W. (3.62) 



/=0 

here 



/ j 22fc-2'(A;-/)! 



P^{X) = ^exp[-(X2-X2)]{exp(-X2)g„(X) + 0FX(l + erf (X))}, (3.63) 



n2n( ]\2 n -i n cy2kh.\ 

- ^I:|---s||-^^ (3-) 



and the normalization constant is equal to 
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Figure 9: Phase distributions for the displaced number state with n = 2 and cto 
Meaning of the hues is the same as in fig. 3. 



= 1 + exp(-X2) I [Q„(X) - 1] d^ 

— 7r 



(3.65) 



The X variable in this case is 



1-s 



ao cos 9, 



(3.66) 



and we have assumed ckq as real. Despite its more complex structure, formula (3.62) 
contains phase distributions Pn{X) that exhibit the main features of the previous phase 
distributions P^^\9), i.e.. cq. (3.3) for a coherent state and eq. (3.36) for a squeezed state. 
Displaced number states have the following Fock expansion: 



\oio,no) = ^6„e^'^"|n), 



(3.67) 



where the amplitudes 6„ and phase-factors 0„ are 

K = (^l^(laol)l'^o) 

1/2 



= exp(-i|aor) (-l)"+-"|aor+-"-L::+-"-(l«or), (3.68) 



n_ = min{n, no}, 

n+ = n + no — ri- = max{n, no}. 



0« 



(n - no) ArgctQ = (n - no)i?o, 



(3.69) 
(3.70) 
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which on insertion into cq. (2.36) give explicitly the Pegg-Barnett distribution P{6). 

Both for coherent states and squeezed states, there was no qualitative difference be- 
tween various phase distributions. Thus, one could say that at least qualitatively, all 
the phase distributions carried the same phase information. Here, we give an example 
of states for which the above statement is no longer true. These are displaced number 
states. The phase properties of such states were discussed earlier by Tanas, Murzakhme- 
tov, Gantsog and Chizhov [1992]. It was shown that there is a qualitative difference 
between the Husimi phase distribution on one side, and the Pegg-Barnett and Wigner 
phase distributions on the other. There is an essential loss of information in the case 
of the Husimi phase distribution. The differences can be interpreted easily when the 
concept of the area of overlap in phase space introduced by Schleich and Wheeler [1987] 
is invoked. Formula (3.62) provides the possibility of deeper insight into the structure 
of the s-parametrized phase distributions. The phase distribution P^^\9) is a result of 
competition between the functions Pn{X), which are peaked at ^ = 0, and the functions 
{Xq — X'^y, which have peaks at ^ = ±7r/2. For s = — 1 only the term with n — k — 
survives, and there is no modulation due to the (—1)"-^ factor. This is the reason why 
the Husimi phase distribution can have at most two peaks, no matter how large is n. 
Both for the Pegg-Barnett phase distribution and P'^'^^O) there are n+l peaks. It is also 
worth noting that despite the fact that the Wigner function yV^^\ [eq. (3.61)] oscillates 
between positive and negative values, the Wigner phase distribution P^°\9) [eq. (3.62)] 
is positive definite. An illustration of the differences between the phase distributions for 
the displaced number states with n = 2 and |«oP = 9 is shown in fig. 9. It is seen that 
the Pegg-Barnett phase distribution is very close to P^^\9), and that they carry basically 
the same phase information, while there is an essential loss of phase information carried 
by P(-i)(^). The Pegg-Barnett and P^^\0) distributions are very similar for given n, 
while P^~^\9) has at most two peaks that become broader as n increases. This example 
shows the difference between a "pure" canonical phase distribution such as the Pegg- 
Barnett distribution, which could be associated with a "pure" phase measurement, and a 
"measured" phase distribution such as P^~^\9), which can be associated with the noisy 
measurement of the phase. The noise introduced by the measurement process reduces the 
phase information that can be inferred from the measured data. 

4 Phase properties of two-mode optical fields 

The single-mode version of the Pegg-Barnett phase formalism can be extended easily into 
the two-mode fields (Barnett and Pegg [1990]) that are often a subject of consideration in 
quantum optics. This leads to the joint phase probability distribution for the phases of the 
two modes, and allows the study of not only the individual mode phase characteristics 
discussed above but also essentially two-mode phase characteristics such as correlation 
between the phases of the two modes. The phase properties of a two-mode field are 
simply constructed from the single-mode phases (see §2.1). The two- mode joint phase 
distribution is given by 




(4.1) 
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This phase distribution can be apphed, similar to the single-mode case, to calculations 
of the mean values of the phase-dependent quantities, such as individual phases, their 
variances, etc. We are often interested not in the individual phases corresponding to either 
mode, but rather in the operators or distributions representing the sum and difference of 
the single-mode phases, which can also be calculated using the joint phase distribution 
[eq. (4.1)]. However, the phase sum and difference values will cover the An range, and 
the integrations over the phase sum and difference variable should be performed over the 
whole range. This approach, although fully justified, is not compatible with the idea 
that the individual phase should be 27r-periodic, and there should be a way to cast the 
phase sum and difference into the 27r range. Such a casting procedure was proposed 
by Barnett and Pegg [1990]. The two approaches, however, give different values for 
the phase sum and difference variances, for example, and one should be aware of the 
differences. Sometimes the original calculations based on the joint phase distribution (4.1) 
have a more transparent interpretation, especially when one considers the intermode phase 
correlations. We shall adduce here examples of both approaches (the quantities obtained 
with the use of the casting procedure will be distinguished by the subscript 27r). The 
casting procedure is described briefly below. 

The possible eigenvalues of the phase-sum operator are 

27r 

where m+ = 0, 1, . . . , 2a, and the eigenvalues of the phase-difference operator are 

27r 

Om. = Ooa-Oob + ^-m_, (4.3) 
cr + 1 

where m_ = —a, —a + 1, . . . ,a. It is seen that the eigenvalue spectra (4.2)-(4.3) of the 
phase sum and difference operators have widths of Att. Since phases differing by 27r are 
physically indistinguishable, the phase sum and difference operators and distributions 
should be cast into a 27r-range (Barnett and Pegg [1990]). The casting procedure can be 
applied to the joint continuous-phase distribution P47r(^+, 0-), defined as 

= l}^{^)'\{(^m^,(^mjf)\', (4.4) 

where 

e± = ea±eb. (4.5) 

As was stressed by Barnett and Pegg [1990], there are many ways to apply the casting 
procedure. However, if the distribution P47r(^+,^_) is sharply peaked, we must avoid 
splitting the original single peak into two parts, one at each end of the 27r-interval. Such 
a poor choice of the 27r-range leads to the same interpretation problems encountered for a 
poor choice of in the single-mode case (Barnett and Pegg [1989]). The casting procedure 
can be applied as follows: 

P2Ae+,e.) = p^^{e+,e.) + p^,{e+ + Si,e. + 62), (4.6) 
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where the shifts 5i and 82 are dependent on the values of ^_ and 



I. 


Si 


= 27r, ^2 = 










for 


0+ 


— \*^U-|- 1 ' 5 U-|- 


+ Stt), 


^_ 




- 6I0+ - 27r| + 6lo_ - 27r, 6lo_ - tt 


11. 




= 0, ^2 = 27r 










for 




e {9o_ - 71, 9o_ 


+ 7r), 


^+ 


e(|^_ 




III. 




= -2n, 62 = 










for 


0+ 


e (^0+ + TT, ^0+ 


+ 37r), 


^_ 




+ 7r,27r- |^+-^o+ -27r|+^o_ 


IV. 




= 0, ^2 = -27r 










for 


9- 


e {9o_ - TT, 


+ 7r), 






+ 37r,47r- -^o_| +^0+)- 



(4.7) 



This analysis of four regions in the (^^+, ^_)-plane to be cut and shifted is close to the 
original idea of Barnett and Pegg [1990], and can be easily understood in a geometri- 
cal representation of the variable transformation. Moreover, as a further consequence 
of the 27r-periodicity of eq. (4.6), one can keep the same shifts 5i and 82 in the whole 
(^+, ^_)-plane without distinguishing any regions. Let us only mention some of the pos- 
sible simplified castings: 

P2.{9+,9_) = P^,{9+,9_) + P,^{9+,9_ ±27r) 
= P,^{9+,9.) + P,^{9+ ±271, 9.) 
= P^^{9+, 9.) + P^^{9+±7r, 9. ±7r) 

(4.8) 

and combinations of the shifts satisfying the condition |5i| + |52| = 27r or | — 1^21 1 = 27r. 
The resulting joint distribution P2-k{9+,9-) is 27r-periodic in 9+ and 9-. Alternatively, 
one can apply the casting procedure to phase distribution (4.1): 



P-2t,{9+,9^) — - 



, (4.9) 



The factor 1/2 occurring in eqs. (4.4) and (4.9) comes from the Jacobian of the transfor- 
mation (4.5) for the variables. The marginal mod(27r) phase-sum, P2-n-{9+), and phase- 
difference, P2n{9-), distributions are given by 

P2.{9^) = / P2.{9+,9_)d9^, (4.10) 

where 

^o± = ^o±±7r. (4.11) 

In the above approach, the casting was prior to the integration. There is another equiv- 
alent manner of obtaining mod(27r) marginal phase sum and difference distributions in 
which the casting is applied after integration. In this approach (Barnett and Pegg [1990]), 
one starts from eq. (4.4) to calculate the 47r-periodic marginal distributions P47r(^±): 

27r-le+-0o+-27r|+eo_ 

Pa.{9+) = / P^^{9+,9_)d9_, (4.12) 

\9+-eo.-2-K\+eo_-2-K 
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47r-|6i--6io_|+6»o+ 

PaAO-) = / p^^{e+,e.)de+. (4.13) 

\e--0o_\+eo+ 

Contrary to the former approach, the casting procedure is now apphed to the single-mode 
distributions P4t:{0±) (Barnett and Pegg [1990]): 

P (ft \ = j P4A0+) + P4Ad+ + 2n) if 9o^+n< 9+ < 00^ + 277 

^2A^+) I p^^(^^) ^ p^^(^^ _ 27r) if +27r<e+< + Stt ^ ' 

and 

p (n ^ i PUO-) + i^4.(^- + 27r) if ^^o_ - TT < ^_ < u ,^ 

~ \ PUO-) + - 27r) if ^^o_ < ^- < + TT . ^ '''' 

Again, due to the 27r-periodicity of P27r(^±) in one can simphfy the recipes (4.14) and 
(4.15) to one of the forms: 

P2.(^±) = P4.(^±)+P4.(^± + 27r) 

= P4.(^±) + P4.(^±-27r) (4.16) 

in the whole intervals 9o± <6±< 9o± + 27r. 

One can analyze analogously the two-mode s-parametrized phase distributions. Here 
we give only one expression for the mod(27r) s-parametrized phase-difference distribution 
for arbitrary density matrix p and any s: 

-1 oo n n 
n=0 A:=0 1=0 

X exp[i{k - l)e^]{l,n - l\p\k,n - k) (4.17) 

with the coefficients G^^^k,!) given by eq. (2.60). Also, by putting G'^^\k,l) 1, 
the mod(27r) Pegg-Barnett phase-difference distribution is obtained as derived by Luis, 
Sanchez-Soto and Tanas [1995]. 



4.1 Two- mode squeezed vacuum 

Single-mode squeezed states, discussed in §3.3, differ essentially from the two-mode 
squeezed states discussed extensively by Caves and Schumaker [1985] and Schumaker 
and Caves [1985]. The Pegg-Barnett phase formalism was applied by Barnett and Pegg 
[1990], and by Gantsog and Tanas [199 Ig] to study the phase properties of the two-mode 
squeezed vacuum, and some of the results are adduced here. 

The two-mode squeezed vacuum state is defined by applying the two-mode squeeze 
operator S{r, if) on the two- mode vacuum, and is given by (Schumaker and Caves [1985]): 

|0,0)(,,^) = S{r,<p)\0,0) 

= (cosh r)'^ exp (e^'^'tanh r aja^) |0,0) 

oo 

= (coshr)-^ ^ (e2'^tanhr)"]n,n), (4.18) 

n=0 
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Figure 10: The joint probability distribution P{9i,92), eq. (4.19), for the two- mode 
squeezed vacuum with r — 0.5. 

where a\ and 02 are the creation operators for the two modes, r (0 < r < 00) is the 
strength of squeezing, and ip (— 7r/2 < 99 < 7r/2) is the phase (note the shift in phase by 
7r/2 with respect to the usual choice of ip). 

The state (4.18), when the procedure described earlier is applied to it, leads to the joint 
probability distribution for the phases 9i and $2 of the two modes in the form (Barnett 
and Pegg [1990]): 

P{ei,92) = (47r2coshV)-^ (1 + tanhV -2tanhrcos(^i + ^2))~^ (4.19) 

One important property of the two-mode squeezed vacuum, which is apparent from 
eq. (4.19), is that P{6i,92) depends on the sum of the two phases only. Integrating 
P{9i, 6*2) over one of the phases gives the marginal phase distribution P{9i) or ^(6*2) for 
the phase 9i or 92'. 

P{9,) = I Pi9,,92)d92 = P{92) = — , (4.20) 

— TT 

meaning that the phases 9i and 92 of the individual modes are distributed uniformly. This 
gives 

TT 

^ 'P+ I 9iP(9i)d9i = (I.,,) = ^, (4.21) 

—IT 

and 

(l>e, + = 2ip, {^0, - = 0. (4.22) 

Thus, the phase-sum operator is related to the phase 2(fi defining the two-mode squeezed 
vacuum state (4.18). 

The two-mode squeezed vacuum has very specific phase properties: the individual 
phases as well as the phase difference are random, and the only non-random phase is the 
phase sum. 
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Figure 11: The phase variances V12 = ([A($e^ + $^2)]^), eq. (4.29), and V+ = ([A($ei + 
^02)V)2it: eq. (4.30), and the phase correlation function C12, eq. (4.28), against the 
squeeze parameter r for the two-mode squeezed vacuum. 

Figure 10 shows an example of the joint phase probability distribution P{6i, ^^2)- The 
ridge, which is parallel to the diagonal of the phase window square reflects the dependence 
of P (^1,^2) on 9i + 92 only. 

The phase distribution P{9i, 62), [eq. (4.19)] is an explicit function of the phase sum, 
but not of the phase difference. This suggests expression of cq. (4.19) in new variables 
{6-^-, 9 ). After applying the casting procedure (see introduction to §4) the joint mod(27r) 
phase distribution is (Barnett and Pegg [1990]): 

P2^(^+,^-) = (47r^coshV)"^ (l+tanhV-2tanhrcose+)~\ (4.23) 
whereas the marginal phase distributions are 

P2^{9+) = (27rcoshV)-^ (l + tanhV-2tanhrcos^+)"\ (4.24) 

P2.{e-) = ^. (4.25) 

The uniform shape of function (4.25) signifies randomness of the phase difference in 
the field [eq. (4.18)]. If the casting procedure is not applied, the marginal distributions 
P{9±) = P^^{9±) have more complicated structures (Barnett and Pegg [1990]). In par- 
ticular, P4ni9-) is not uniform because of the integration limits in eq. (4.13). In general, 
the mod(47r) distribution has no unique shape signifying randomness of the phase sum 
or difference. There are many distributions in the 47r-range leading to a flat mod(27r) 
function. 

The two-mode variance of the phase-sum operator can be calculated according to the 
general formula: 

{[A{^e, + ^9,)?) = ((A|.,j2) + ((A|.,j2)+2Ci2 (4.26) 
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in terms of the individual phase variances ((A^^i^^)^) the phase correlation function 
(correlation coefficient) 

Ci2 = ($61^62) - (^ei) (^62) 

TT 77 TT TT 

= 11 e,e2P{ei, 62) dOi - ( / ^iP(^i) d^i) ( J ^2P(^2) d^2) ■ (4.27) 

— TT — TT — TT — TT 

The variances ((A$6i^ 2)^) are simply 7r^/3 [because of eq. (4.20)], and the phase correlation 
function C12 is equal to: 



(tanl 



C12 = -2(cosh r)-^ y ^7'"' = -2dilog(l-tanhr). (4.28) 

^fc (ri-A;)2 

This correlation function describes the correlation between the phases of the two modes of 
the two- mode squeezed vacuum. In fig. 11 the correlation coefficient as well as the phase 
variances are plotted against the squeeze parameter r. The correlation is negative and, as 
r tends to infinity, approaches — tt^/S asymptotically. Finally, phase variance (4.26) has 
the form: 

{[A{^e^+K)?) = 2^-4dilog(l-tanhr). (4.29) 

The strong negative correlation between the two phases lowers the variance (4.29) of the 
phase-sum operator. For r — 00, this variance tends asymptotically to zero, which means 
that for very high squeezing the sum of the two phases becomes well defined (phase locking 
effect) . 

The ( "single- mode" ) mod(27r) phase-sum variance is (Barnett and Pegg [1990]): 



i+k 



n 



— TT 
,2 



= — + 4dilog(l+tanhr). (4.30) 
3 

As the squeezing parameter r varies from to 00, the mod(27r) variance [eq. (4.30)] de- 
creases from 71^/3 to zero, whereas the two-mode phase-sum variance [eq. (4.29)] changes 
from 27rV3 to zero with increasing r. Hence, both variances (4.29) and (4.30), reveal the 
fact that the phase sum becomes perfectly locked in the limit of large squeezing (r 00). 
The value of the variance (4.30) describes random phase sum for zero squeezing. 
In this case of r = 0, the two-mode variance (4.29) is twice as much as the mod(27r) 
phase-sum variance (4.30), since it shows randomization of the two phases, and ^$2^ 
separately. 

As was stressed in §4, both the original distributions, given by eqs. (4.19) and (4.20), 
and the mod(27r) distributions, given by eqs. (4.23)-(4.25), are valid and useful. However, 
some care is required when interpreting the results obtained in both ways. The phase- 
sum variance has generally different values, as seen from fig. 11, in the two approaches. 
The original distributions are better for understanding the intermode phase correlation, 
which can be calculated explicitly from eq. (4.27), while for the mod(27r) distribution 
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the correlation is concealed in the value of the phase variance (4.30) and is not seen 
explicitly. On the other hand, the mod(27r) results have a clear interpretation for the sum 
and difference of the individual phases treated as single phase variables. 

Generalizing formula (2.15) and taking into account the fact that the two-mode 
squeezed vacuum is a "physical state", we can calculate the expectation values of the 
phase exponential operators in the following way (Gantsog and Tanas [1991g]): 

(exp(imil>0jexp(im2l>e2)) = (e5cp(imi$isG) exp(im2$2SG)) 

oo oo l-\-m 

— (coshr)"^ ^ ^ (^e^'^tanhr) {l,l\n, k){n + rrii, k + m2\m,m) 

n,k=0 rn.l=0 

e'"^tanhr) 5^,,^,, (4.31) 
where for brevity we denote ((...)) = (r,<^)(0, 0|(...)|0, 0)(r,<^). The operators 

oo 

Scp(imi,2$i,2SG) = '^\n){n + mi^2\ (4.32) 

n=0 

are the Susskind-Glogower phase operators for the two modes. Formula (4.31) is strikingly 
simple, and shows that only exponentials of the phase sum have nonzero expectation 
values. 

Using eq. (4.31), the following results for the cosine and sine of the phase-sum operator 
are obtained (Gantsog and Tanas [1991g]): 

(cos(<l>e)^ -I- $5)2)) = tanhrcos2v9, 

(sin(l>ei + I'e2)) = tanhrsin2(^, (4.33) 



(cos2(l>0, + $^2)) 



1 1 

— I — (tanhr)^ cos4(^, 

2 2 



(sin^(l>ei +^62)) = \- ^(tanhr)^cos4(^. 



(4.34) 



([Acos($,,+$,J]2) = ([Asin(|.,,+$,J]2) = l(coshr)-2. (4.35) 

For very large squeezing (r 00, tanhr — )■ 1, coshr 00), the expectation values (4.33) 
and (4.34) of the functions of the phase-sum operator become asymptotically correspond- 
ing functions of the phase 2(f, confirming the relation between the phase sum and 2(fi that 
is already apparent from eq. (4.22). It is interesting that the expectation value of the 
phase-sum operator is equal to 2ip irrespective of the value of r, whereas for the sine and 
cosine functions correspondence is obtained only asymptotically. The variances (4.35) 
then become zero and the sine and cosine of the phase sum are well defined. 

It should, however, be emphasized that the expectation values calculated according to 
the Pegg-Barnett formalism depend on the choice of the particular window of the phase 
eigenvalues. If a choice different from that made above were made, the clear picture of 
the phase properties of the two-mode squeezed vacuum would be disturbed. For exam- 
ple, the value of the correlation coefficient (4.28) would be different, and the phase-sum 
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variance (4.26) would not tend asymptotically to zero. However, formulas (4.31)-(4.35), 
because of the way they have been calculated do not, in fact, depend on the choice of the 
phase window. This gives us the opportunity to make a choice which introduces consis- 
tency in the behavior of the phase itself and its sine and cosine functions. Another way 
of making the choice is to minimize the variance (4.26) of the phase-sum operator. 

4.2 Pair coherent states 

Pair coherent states introduced by Agarwal [1986, 1988] are quantum states of the two- 
mode electromagnetic field, which are simultaneous eigenstates of the pair annihilation 
operator and the difference in the number operators of the two modes of the field. Agarwal 
[1988] has discussed the nonclassical properties of such states, showing that they exhibit 
remarkable quantum features such as sub-Poissonian statistics, correlations in the number 
fluctuations, squeezing, and violations of the Cauchy-Schwarz inequalities. He has also 
presented results for fluctuations in the phase of the field using the Susskind-Glogower 
phase formalism. The phase properties of such states on the basis of the Pegg-Barnett 
formalism were studied by Gantsog and Tanas [1991e], and by Gou [1993]. Phase distri- 
butions for squeezed pair coherent states were analyzed by Gerry [1995]. 

The pair coherent states are defined by Agarwal [1988] as eigenstates of the pair- 
annihilation operator 

&b\C,q) = C\Cq), (4.36) 

where C is a complex eigenvalue and q is the degeneracy parameter, which can be fixed 
by the requirement that |C,?) is an eigenstate of the difference of the number operators 
for the two modes: 

(a^a-Pb)\(:,q) = q\C,q). (4.37) 

The solution to the above eigenvalue problem, assuming q to be positive, is given 
by (Agarwal [1988]): 

oo An 

IC«) = 'V,Ep(;^^l" + ?."). (4.38) 

where Ng is the normalization constant 

SsFT^j = l('KI)-'^.PiKI)l-"^- (4.39) 

The state \n + q,n) is a Fock state with n + q photons in mode a and n photons in mode 
b. If the complex number ^ is written in the form: 

C=|C|exp(i(^), (4.40) 

the state (4.38) can be written as 

oo 

\C,q) = T.bne'''''\n + q,n), (4.41) 

n=0 
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where 

Now, the phase properties of the state (4.41) can be studied easily using the Pegg-Barnett 
formahsm in a standard way as described above. The resulting joint probability distribu- 
tion for the phases 9a and 9^ of the two modes is given by (Gantsog and Tanas [1991e]): 

P{0a,9,) = ^S^l + 2Y^KbkCos[{n-k){9a + 9,)]y (4.43) 



where 6„ is given by eq. (4.42). For q = 0, formula (4.43) can be written in the following 
simple form: 

P(9a,9b) = ^exp[2|C|cos(^„ + ^6)]. (4.44) 

As in the case of the two- mode squeezed vacuum, the joint phase probability distribu- 
tion depends on the sum of the two phases only, which means strong correlation between 
the two phases. Again, the only non-uniformly distributed phase quantity is the phase 
sum 9a + 9b. This suggests re-expression of the phase distribution (4.43) in new variables 
of the phase sum, 6^ = 6a + 9b, and phase difference, 9^ = 9a — 9h- After applying the 
casting procedure, the mod(27r) Pegg-Barnett distribution P27r(^+,^-) takes the form: 

P2.{.9+,9_) = ^|l + 2f;6Acos[(n-A;)^+]l, (4.45) 

and the marginal distributions are: 

P2^{9+) = ^|l + 2f;6Acos[(n-fc)^+]i, (4.46) 

2^ ( n>k J 

and 

P2.{d-) = ^- (4.47) 

For completeness of our discussion and by analogy with our presentation of the single- 
mode models, we now give expressions for various s-parametrized phase distributions. 
Thus, the mod(27r) two-mode s-parametrized phase distribution is equal to 

Pt\0+,0-) = ^ll + 2Y,bnbkG^'\n,k)G('\n + q,k + q)cos[in-k)9+]\, (4.48) 

where the coefficients G^^\n,k) are given by eqs. (2.60)-(2.62). The mod(27r) marginal 
s-parametrized phase-sum distribution is 

Pi^{0+) = —ll + 2Y,bnbkG'^'\n,k)G'^'\n + q,k + q)cos[{n-k)9+]\.{4A9) 
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The mod(27r) s-paramctrized phase-difference distribution P2^'{9^) and the single-mode 
ones, P(*)(6'(j) and P'^^\9b), are uniform: 

pi:V_) = p^'\ea) = p^'\e,) = ^. (4.50) 

The distributions (4.48)-(4.50), similar to the distributions (4.45)-(4.47), reveal the fun- 
damental phase properties of pair coherent states. 

The correlation coefficient Cab, eq. (4.27), [subscripts 1, 2 should be replaced by a and 
b, respectively] is given in this case by the formula: 

where 6„ is given by eq. (4.42). This correlation is negative and lowers the variance of 
the phase-sum operator. For |^| oo, this coefficient approaches — 7r^/3, the phase-sum 
variance becomes zero, and we have the classical situation of perfectly defined phase sum 
(the phase locking effect). This phase correlation coefficient can be contrasted with the 
photon number correlation coefficient, considered by Agarwal [1988], which increases as 
\(\ increases. The sine and cosine functions of the phase-sum operator were also obtained 
by Gantsog and Tanas [1991e] and compared to their counterparts obtained by Agarwal 
[1988], who used the Susskind-Glogower approach. 

4.3 Elliptically polarized light propagating in a nonlinear Kerr 
medium 

To describe propagation of elliptically polarized light in a nonlinear Kerr medium, a 
two-mode description of the field is needed. The quantum nature of the field results 
in the appearance of such quantum effects as photon antibunching (Ritze and Bandilla 
[1979], Tanas and Kielich [1979], Ritze [1980]) and squeezing (Tanas and Kielich [1983, 
1984]). Tanas and Kielich have shown that as much as 98 percent of squeezing can be 
obtained when intense light propagates in a nonlinear Kerr medium. They referred to this 
effect as self-squeezing. Agarwal and Puri [1989] re-examined the problem of propagation 
of elliptically polarized light through a Kerr medium, considering not only the Heisenberg 
equations of motion for the field operators, but also the evolution of the states themselves. 
Quantum fluctuations in the Stokes parameters of light propagating in a Kerr medium 
were discussed by Tanas and Kielich [1990], and by Tanas and Gantsog [1992b]. 

The following effective interaction Hamiltonian can be used to describe the propagation 
of eUiptically polarized hght in a Kerr medium (Tanas and Kielich [1983, 1984]): 

Hi — -Hk fal^Oi -l- 02^02 + 4cia| 020102) , (4-52) 

where oi and 02 are the annihilation operators for the circularly right- ("1") and left- 
("2") polarized modes of the field, k is the coupling constant, which is real and related to 
the nonlinear susceptibility tensor x*^^^ of the medium, and d is the asymmetry parameter 
describing the coupling between the two modes. For a fully symmetrical susceptibility 
tensor, d=\. Otherwise, (i 7^ and describes the asymmetry of the nonlinear properties 
of the medium (Ritze [1980], Tanas and Kiehch [1983, 1984] ). 
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Using the Hamiltonian (4.52), one can obtain the cvohition operator U{t), and as- 
suming that the initial state of the field is a coherent state of the elliptically polarized 
light, one obtains for the resulting state of the field (Agarwal and Puri [1989]): 

IV'(t)) = U{T)\a,,a2) 

= m Kibn2 exp< i(ni(/7i + n2</?2) 

+ i^ [ni{ni - 1) +^2(712 - 1) +4dnin2]||ni,n2), (4.53) 

where r = n{uj)kz/c {n{uj) with the refractive index), and the coefficients 6„^_2 given 
by eq. (3.8) with |q;iP and |Q;2p as the mean numbers of photons for the circularly right- 
and left-polarized modes, respectively, whereas (fi^2 are the phases of the coherent states 
of the two modes. 

The state (4.53) is the two-mode state of the field, and the two-mode generalization 
of the Pegg-Barnett formalism used by Gantsog and Tanas [1991c] leads to the following 
joint probability distribution for the continuous-phase variables, 9i and $2, of the two 
modes: 



P{ei,e2 



m m bnibn2 expj -ini6i - m202 

1 r\m n 



ni=0 712=0 



+ i^[ni(ni - 1) + ^2(712 - 1) + 4cinin2]| 



(4.54) 



The phase distribution function P{9i,92) describes the phase properties of elliptically 
polarized light propagating through a Kerr medium, which were discussed in detail 
by Gantsog and Tanas [1991c]. Figure 12 shows an example of the evolution of P{9i, 62). 
It is seen that the peak is shifted and broadened during the evolution. Since the numbers 
of photons in the two modes are different, one can see that the shift of the peak and its 
broadening is asymmetric. The intensity-dependent phase shift is bigger for the mode 
with higher number of photons. This corresponds to the classical effect of self-phase mod- 
ulation in a nonlinear Kerr medium. The quantum description shows not only the shift 
but also the broadening of the phase distribution (phase diffusion). 

Integration of the distribution function P{6i,62) over one of the phases 61 or 62 leads 
to the marginal distribution ^(6*2) or ^(6*1) for the individual phases. All single-mode 
phase characteristics of the field can be calculated using these distributions, and the 
corresponding formulas were given by Gantsog and Tanas [1991c]. 

In addition to the phase properties of the individual modes, it is interesting, in the 
two-mode case, to study the behavior of the phase difference between the two modes. In 
the Pegg-Barnett formalism, the phase-difference operator is simply the difference of the 
phase operators for the two modes, so the mean value of the phase-difference operator is 
the difference of the mean values of the single-mode phase operators. To calculate the 
variance of the phase-difference operator, we can use the relation: 

{[A{^e, - ^02)?) = ((A|.,J^) + ((A|.,J^) - 2Ci2, (4.55) 

where the last term is the correlation coefficient between the phases of the two modes and 
can be calculated by integration of P{9i,92) according to eq. (4.27). Thus, the resulting 
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Figure 12: Evolution of the joint probability distribution P{6i,62), eq. (4.54), of light 
propagating in a Kerr medium: |q!iP = 0.25, |q;2|^ = 4 and o? = 1; (a) r = 0, (b) r = 0.1, 
(c) T = 0.2, (d) r = 0.3. 

formula is (Gantsog and Tanas [1991c]): 

CMr) = E E /12/21 - f E E hibl] ( E E h^bl) , (4.56) 

ni>nj n2>n2 \ni=n^ n2>n2 / \n2=n2ni>n^ / 

where 

/i,- = 2bnXr — ^^sin<^-(ni-n^) + - 1 + 2d(n,- + n;) [. (4.57) 

A graphic illustration of the correlation function (4.56) is shown in the left-hand panel of 
fig. 13. The strength of the correlation depends crucially on the value of the asymmetry 
parameter d. The highest values of the correlation are obtained for d — 1/2. This means 
that the minimum of the phase-difference variance, in view of eq. (4.55), is obtained for 
d = 1/2. The phase- difference variance is shown in the right-hand panel of fig. 13. It 
was shown (Tanas and Gantsog [1991]) that, similar to the single-mode case, dissipation 
destroys the periodicity of the evolution and broadens the phase distribution. 

Recently, the phase properties of light propagating in a Kerr medium have been re- 
considered (Luis, Sanchez-Soto and Tanas [1995]) from the point of view of the Hermitian 
phase-difference operator introduced by Luis and Sanchez-Soto [1993b, 1994], which is 
based on the polar decomposition of the Stokes operators. This example shows clearly the 
difference between the Pegg-Barnett and Luis-Sanchez-Soto phase-difference formalism, 
which is most visible for weak fields. The Luis-Sanchez-Soto phase-difference operator 
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Figure 13: Evolution of the intermodc phase correlation function Cyi{t\ eq. (4.56), 
and the phase-difference variance ([A(<l>e^ — $02)]^)) ^q. (4.55), of light propagating in a 
Kerr medium. Thin solid line: — jccip = 0.25, \ol2^ — 4 and d — \\ bold solid line — 



0.25, |a2| = 4 and d 
d=\\ bold dashed line — |q;i| 



1/2; thin dashed hnc — \<^\^ = 0.25, |Q;2p = 0.25 and 
= 0.25, |a2p = 0.25 and d = 1/2. 



differs from the Pegg-Barnett phase-difference operator, which is simply the difference of 
the phase operators of the two modes. For strong fields both formalisms give the same 
results. The nonlinear Kerr medium appears to be a good testing ground for different 
phase approaches. 

As shown by Gantsog and Tanas [1991a], superpositions with any number of compo- 
nents can be obtained in the process of light propagation in the Kerr medium (similar 
to the anharmonic oscillator model described in §3.5) if the evolution time r is taken 
as a fraction M/N of the period, where M and N are mutually prime integers. Exact 
analytical formulas for finding the superposition coefficients were given for any N. The 
joint phase probability distribution P{9i, 62) splits into separate peaks if the state of the 
field becomes a discrete superposition of coherent states, and this is a very spectacular 
way of presenting such superpositions. Some examples are shown in fig. 14. 

4.4 Second-harmonic generation 

Second-harmonic generation is probably the best known nonlinear optical process. In the 
quantum picture we deal with a nonlinear process in which two photons are annihilated 
and single photon with doubled frequency is created. The quantum states of the field 
generated in the process exhibit a number of unique quantum features such as photon 
antibunching (Kozicrowski and Tanas [1977]) and squeezing (Mandel [1982], Wu, Kim- 
ble, Hall and Wu [1986]) for both the fundamental and second- harmonic modes (for a 
review and literature see Kielich, Kozierowski and Tanas [1985]). Nikitin and Masalov 
[1991] discussed the properties of the quantum state of the fundamental mode, calcu- 
lating numerically the quasiprobability distribution function Q{a, a*) for it. They sug- 
gested that the quantum state of the fundamental mode evolves, in the course of the 
second-harmonic generation, into a superposition of two macroscopically distinguishable 
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Figure 14: The joint probability distribution P{6i, 62), eq. (4.54), of light propagating in 
a Kerr medium. |q;iP = |Q;2p = 4, r = 27r/2, and (a) d = 0; (h) d = 1/2. 

states, similar to the superpositions obtained for the anharmonic oscillator model (Yurke 
and Stoler [1986], Tombesi and Mecozzi [1987], Miranowicz, Tanas and Kielich [1990], 
Gantsog and Tanas [1991f]), or a Kerr medium (Agarwal and Puri [1989], Gantsog and 
Tanas [1991a]). Gantsog, Tanas and Zawodny [1991a] discussed the phase properties of 
the field produced in the second-harmonic generation process. 

To describe second-harmonic generation, the following model Hamiltonian is used: 

H ^ Ho + Hi = hua^a + 2hujPb + hg{b'^a^ + bd^'^), (4.58) 

where d (a^) and b (6^) are the annihilation (creation) operators of the fundamental mode 
of frequency cu and the second-harmonic mode at frequency 2a;, respectively. The coupling 
constant g, which is real, describes the coupling between the two modes. Since Hq and Hi 
commute, there are two constants of motion: Hq and Hi, Hq determines the total energy 
stored in both modes, which is conserved by the interaction Hi. The free evolution can 
be thus factored out, and the resulting state of the field can be written as 

|*(i)) = exp(-i^iV^)|*(0)), (4.59) 

where 1^1/(0)) is the initial state of the field. Since the interaction Hamiltonian Hi is not 
diagonal in the number-state basis, the numerical method of diagonalization of Hi may 
be applied to find the state evolution (Walls and Barakat [1970]). 

Let us assume that initially there are n photons in the fundamental mode and no 
photons in the second-harmonic mode; i.e., the initial state of the field is |n, 0) = |n)|0). 
Since Hq is a constant of motion, we have the relation: 

(a^a) + 2(6^6) = constant = n, (4.60) 

which implies that the creation of k photons of the second-harmonic mode requires anni- 
hilation of 2k photons of the fundamental mode. Thus, for given n, we can introduce the 
states 

iV'i"^) = \n-2k,k), A;-0,l,...,[n/2], (4.61) 
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Figure 15: Evolution of the joint phase probabihty distribution P{6a,6b), eq. (4.71), in 
the second-harmonic generation. The initial mean number of photons of the fundamental 
beam is |q;oP = 4, and gt is the dimensionless scaled time. 
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where [n/2] denotes the integer part of n/2, which form a complete basis of states of the 
field for given n. We have 

= Kn'kk', (4.62) 

meaning that the constant of motion Hq splits the field space into orthogonal subspaces, 
which for given n have the number of components equal to [n/2] + 1. In such a basis, the 
interaction Hamiltonian has the following nonzero matrix elements: 

= ng^{k + l)(n - 2k){n -2k- 1), (4.63) 

which form a symmetric matrix of the dimension ([n/2] + 1) x ([n/2] + 1) with real nonzero 
elements (we have assumed g real) located on the two diagonals immediately above and 
below the principal diagonal. Such a matrix can be easily diagonalized numerically (Walls 
and Barakat [1970]). 

To find the state evolution, we need the matrix elements of the evolution operator: 

dn,k(t) = (4"^|exp(-i^it/;i)|4")). (4.64) 

If the matrix U is the unitary matrix that diagonalizes the interaction Hamiltonian matrix 
given by eq. (4.63), i.e.: 

U-^Ht'^ij = %xdiag(Ao,Ai,...,A[„/2]), (4.65) 

then the coefficients dn,k{t) can be written as 

[n/2] 

dn,k(t) = 5^ e-'^*^*C/fc,C/o* , (4.66) 

where Aj are the eigenvalues of the interaction Hamiltonian in units of kg. Of course, 
the matrix U as well as the eigenvalues Aj are defined for given n and should have the 
additional index n, which we have omitted to shorten the notation. Moreover, for real 
g the interaction Hamiltonian matrix is real, and the transformation matrix [/ is a real 
orthogonal matrix, so the asterisk can also be dropped. 

The numerical diagonalization procedure gives the eigenvalues A, as well as the ele- 
ments of the matrix U, and thus the coefficients d„j,.{f?) can be calculated according to 
eq. (4.66). It is worthwhile to note, however, that due to the symmetry of the Hamil- 
tonian the eigenvalues Aj are distributed symmetrically with respect to zero, with single 
eigenvalue equal to zero if there is an odd number of them. When the eigenvalues are 
numbered from the lowest to the highest value, there is an additional symmetry relation: 

UkiUoi = {-'^TUk,[n/2\-iUo,[n/2\-i , (4.67) 

which makes the coefficients dn,k{t) either real {k even) or imaginary {k odd). This prop- 
erty of the coefficients dn^kit) is very important and in some cases allows exact analytical 
results to be obtained. 
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With the coefficients (i„ ^(t) available, the resulting state of the field (4.59) can be 
written, for the initial state |n, 0), as 

\^^-\t)) = E^fcWI^i"^). (4.68) 

k=0 

The typical initial conditions for the second-harmonic generation are: a coherent state of 
the fundamental mode and the vacuum of the second-harmonic mode. The initial state 
of the field can thus be written as 

oo 

IV'(O)) = J2^n\n,0), (4.69) 

n=0 

where c„ = 6„e'"''*''' is the Poissonian weighting factor (3.8) of the coherent state |ao) with 
the phase (pa — ArgctQ. With these initial conditions, the resulting state (4.59) is given 

by 

oo oo [n/2] 

m)) = Y.(^n\^^''\t)) = Y.CnJ2dnAt)\n-2k,k). (4.70) 

n=0 n=0 fc=0 

Equation (4.70), describing the evolution of the system is the starting point for a further 
discussion of second-harmonic generation. If the initial state of the fundamental mode is 
not a coherent state, but has a decomposition into number states of the form (4.69) with 
different amphtudes Cn, eq. (4.70) is still valid if appropriate taken. This is true, 

for example, for an initially squeezed state of the fundamental mode. The coefficients 
dn,k{t) have been calculated numerically to find the evolution of the field state (4.59), and 
consequently, its phase properties (Tanas, Gantsog and Zawodny [1991a, 1991b], Gantsog, 
Tanas and Zawodny [1991a]). 

Repeating the standard procedure of the Pcgg-Barnctt formalism with the field 
state (4.59), the joint phase probability distribution is obtained in the form: 



(27r) 



oo [n/2] 

dn,k{t) exp{-i[(n - 2k)9a + k9b - k{2ipa - (pb)]} 

n=0 k=0 



(4.71) 



where 6a and 6h are continuous-phase variables for the fundamental and second-harmonic 
modes, and the phases (fa and iff, are the initial phases with respect to which the dis- 
tribution is symmetrized. It is interesting, that formula (4.71) depends, in fact, on the 
difference 2(pa — <Pb: which reproduces the classical phase relation for second-harmonic gen- 
eration. Classically, for the initial conditions chosen here, this phase difference takes the 
value n/2, which turns out to be a good choice to fix the phase windows in the quantum 
description as well. 

The evolution of the joint probability distribution P{9a, 0^), given by eq. (4.71), is il- 
lustrated graphically in fig. 15. At the initial stage of the evolution the phase distribution 
in the 9a direction (fundamental mode) is broadened, while a peak of the second-harmonic 
mode phase starts to grow. The emergence of the peak at = confirms the classical 
relation 2(pa — (fb = 7r/2, which has been applied to fix the phase window. The phase 
distribution in the 9^ direction narrows at the beginning of the evolution, meaning less 
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uncertainty in the phase of the second harmonic. However, for later times the distribu- 
tion P{Oa, Ob) sphts into two peaks, which resembles the splitting of the a*) function 
found by Nikitin and Masalov [1991]. For still later times, more and more peaks appear in 
the distribution P{9a,0h), and this distribution becomes more and more uniform, which 
means randomization of the phase. The route to the random phase distribution, how- 
ever, goes through a sequence of increasing numbers of peaks. The splitting of the joint 
phase distribution can be understood if one realizes that the mean number of photons 
of the second harmonic oscillates and after reaching the maximum the second-harmonic 
generation becomes, as a matter of fact, the down-conversion process which exhibits a 
two-peak structure of the phase distribution in the direction of the fundamental mode 
(see §4.5). The appearance of new peaks may thus be interpreted as a transition of the 
process from the second-harmonic to the down- conversion regime, and vice versa. The 
phase variances for both modes tend asymptotically to the value 7r^/3 of the randomly 
distributed phase (Gantsog, Tanas and Zawodny [1991a]); however, it has turned out that 
partial revivals of the phase structure can be observed during the evolution (Drobny and 
Jex [1992]). 

It is also interesting to study the phase distribution of the field produced by second- 
harmonic generation with other than coherent initial states of the fundamental mode. 
Such studies were performed by Tanas, Gantsog and Zawodny [1991b], showing for exam- 
ple that even for a second harmonic generated by an initial number state the joint phase 
probability distribution has a modulation structure owing to the intermode correlation 
that develops in the process of the evolution. 

4.5 Parametric down conversion with quantum pump 

The parametric down-conversion process with quantum pump, which is a subharmonic 
generation process, can be described by the same model Hamiltonian, [eq. (4.58)] as 
the second-harmonic generation. The initial conditions distinguishing the two processes 
are the following: For the subharmonic generation process mode b is initially populated 
while mode a is in the vacuum state. The distinction between the two processes is far 
from trivial, and the states generated in the two processes have quite different properties 
(Gantsog, Tanas and Zawodny [1991b], Jex, Drobny and Matsuoka [1992], Tanas and 
Gantsog [1992a, 1992b], Gantsog, Tanas and Zawodny [1993]). 

Let us assume, in analogy to our analysis of second-harmonic generation, that initially 
there are n photons in the pump mode (6) and no photons in the signal mode (a); i.e., 
the initial state of the field is |0, n) — \0)a\n)b. Since Hq is a constant of motion, we have 
the relation 

{a^d) + 2{b^b) = constant = 2n, (4.72) 

which implies that the annihilation of k photons of the pump mode requires creation of 
2k photons of the signal mode. Thus, for given n, we can introduce the states 

IV'S) = \2k,n-k), fc = 0,l,...,n, (4.73) 

which should be compared to the corresponding expression (4.61) for the second-harmonic 
generation. 
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Figure 16: Evolution of the joint phase probabihty distribution P{9a,0i,), eq. (4.76), in 
the process of parametric down conversion with quantum pump. The initial mean number 
of photons in the mode 6 is |/9oP = 4. 
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Figure 17: Evolution of the phase variances {{A^gJ"^), eq. (4.77), and ((A$5iJ^), eq. 
(4.78), in the parametric down conversion. The initial mean number of photons in the 
mode b is |/3oP = 4. 



Proceeding along the same lines as in second-harmonic generation, the resulting state 
of the field can be written as 

oo n 

m)) = E^ne'"^''Ec^2n,fe(i)|2A;,n-A;), (4.74) 

n=0 k=0 

where the coefficients d2n,kif) are given by 

d2n,k{t) = (2A;,n- A;|exp(-i^ii) |0,n), (4.75) 

whereas now the c„ = fc^e"*'^'' are the Poissonian weighting factors (3.8) for the initially 
coherent state |^o = \^q\ exp (if/^b)) of the mode h. Again, the method of numerical diago- 
nalization is used to calculate the coefficients d2n,k{t) and, in effect, the phase properties 
of the state (4.74). 

The joint phase probability distribution in this case is given by 



(27r)^ 



^KY^ d2n,k{t) exp{-i[2/c6'a + (n - k)9b + k{2ipa - (pb)]] 

n=0 ifc=0 



. (4.76) 



As for the second-harmonic generation, we similarly take 2ipa — </96 = 7r/2 to fix the phase 
windows. The evolution of the joint probability distribution P(^a, 6b) for parametric down 
conversion with the mean number of photons |/3oP = 4 is shown in fig. 16. Comparison 
of figs. 15 and 16 shows immediately a striking difference between the phase properties 
of the fields obtained in the two processes. The state produced in the down-conversion 
process acquires from the very beginning the two-peak structure in the 6a direction, which 
suggests the appearance of a superposition of two states in the resulting field. The two 
peaks which appear at the beginning of the evolution correspond, in fact, to the two-peak 
phase distribution of the squeezed states (see §3.3). At later stages of the evolution ran- 
domization of the two phases takes place, similarly as for the second-harmonic. However, 
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Figure 18: Phase distribution P{9a), eq. (4.79), for gt — 0.3, in the parametric down 
conversion. 



the symmetry with respect to 6a is preserved. The two-peak structure of the phase dis- 
tribution has already appeared, although not in its pure form, in the phase distribution 
for second harmonic generation (fig. 15). Its appearance can be ascribed to the down- 
conversion process that has overcome second-harmonic generation at this stage of the 
evolution. The transition from the one-peak phase distribution to the two-peak distribu- 
tion makes a qualitative difference between the two field states, and is a sort of "phase 
transition" . 

Once the joint phase distribution P(^a, ^fc) is known, all quantum-mechanical phase 
expectation values can be calculated. In particular, the phase variance for the signal mode 
can be calculated according to the formula: 



TT 



Re ^ hnhn' 



n>n 



cxp[-i(n - nO(2(/?a - </?&)] 
{n — n')'^ 



k=Q 

and for the pump mode we have 

TT TT 
— TT — TT 



(4.77) 



+ 4Re J2 bnbn' _ ^/)2 d2n,k{^)d*2^',kit)^ 



n>n' 



(4.78) 



k=0 



where we have used eq. (4.76), and we take 2(pa — (fb = 7r/2. The time evolution of the 
phase variances can be calculated numerically using these expressions for given initial 
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field states. The dynamical behavior of the phase variances calculated from eqs. (4.77) 
and (4.78) is illustrated graphically in fig. 17 for |/3oP = 4. The dashed line 7r^/3 marks 
the variance for the state with random distribution of phase. It is apparent that the phase 
variance of the signal mode starts from the value 7r^/3, dips into the minimum, and after 
a few oscillations again becomes close to 7r^/3. For comparison, the phase variance for 
the ideal squeezed state is also shown. The two variances are initially indistinguishable, 
but the phase variance for the squeezed state approaches monotonically its asymptotical 
value 7r^/4, while for the quantum- pump case the phase variance of the signal mode begins 
to oscillate at later times. This confirms the statement that there a limit is imposed by 
the quantum fluctuations of the pump on the applicability of the ideal down-converter 
model. The phase variance of the pump mode increases rapidly from its initial value 
for the coherent state, and also shows oscillatory behavior approaching the value 7r^/3 
at the long-time limit. Thus, the long-time effect of the quantum fluctuations of the 
pump mode is the randomization of the phase distribution for both signal and pump 
modes. This randomization process is not monotonous, and it turned out that at least 
partial revivals of the phase structure are possible during the evolution (Gantsog [1992], 
Gantsog, Tanas and Zawodny [1993]). 

Integrating P{Oa, Of,) over one of the phases leads to the marginal phase distributions 
P{9a) and P{9b) for the phases 9a and 9b of the individual modes. We have 

-I ( n n' 

m) = ^ l + 2ReEMn'EE^2n,fe(t)c?L',/c'W 
I n>n' k=Ok'=0 

X exp [-i{k - k'){29a + 2(pa - iPb)] Sn-n',k-k' \, (4.79) 



PiOb) = — l + 2ReE&n6n'EVfeW^2<feW 

I n>n' k=0 

X exp[-i{n-n)9b]j- (4.80) 

The phase distribution P{9a) for the signal mode is shown in fig. 18 for gt = 0.3; i.e., for 
the time at which the squeezing in the signal mode has its maximum value. For comparison 
we show the phase distribution for the squeezed vacuum for r = 2\/3Q\gt = 1.2. The effect 
of quantum fluctuations of the pump is seen as the broadening of the phase distribution 
with respect to that for the ideal squeezed state. 



5 Conclusion 

In this article we have reviewed some recent results concerning the quantum-phase de- 
scription of optical fields. We have focused our attention on the real fields that can be 
generated in practice in various nonlinear optical processes. So, we rather avoided dis- 
cussions of the phase formalisms as such and tried to exploit their practical applicability 
in the description of optical fields. In the description of the phase properties we used two 
different, though related formalisms: the Pegg-Barnett Hermitian phase formalism and 
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the formalism based on s-parametrized phase distributions. The Pegg-Barnett Hermitian 
phase formahsm is a good example of the concept of the phase as a physical property of a 
single field mode represented by a Hermitian phase operator canonically conjugate one to 
the number operator. It allows one to obtain the phase distributions for the fields, mean 
values and variances of the phase, and other phase characteristics of the field in a rea- 
sonably simple way, both from the conceptual as well as the calculational point of view. 
The phase distributions obtained from this formalism are 27r-periodic, positive definite 
and normalized. They can be treated as a good representation of the quantum state of 
the field and can be referred to as canonical phase distributions. 

Another description of the optical phase used by us is that based on the s-parametrized 
quasiprobability distributions, which can give phase distributions that can be both nar- 
rower and broader (depending on s) than the Pegg-Barnett phase distribution, but these 
distributions with s < — 1 can be associated with some noisy, real measurements of the 
phase probability distribution and can be referred to as measured phase distributions. Us- 
ing the examples of real-field states presented here, we tried to show the similarities and 
differences that one encounters when various phase distributions are applied to describe 
a particular field state. Our choice of the field states is, of course, a bit arbitrary, and we 
relied to a large extent on our own results. We believe, however, that our review covers a 
number of field states important for quantum optics, and that the results presented here 
may prove interesting. We have also attempted to give a more or less complete review of 
the literature on the subject, but the subject of quantum phase is still a "hot" one and 
the literature is growing rapidly. 
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Appendix A. Garrison- Wong phase formalism 

Garrison and Wong [1970] constructed the phase operator I»gw using the relation: 



for any g, f & "H^, where 'H? is the Hilbert space in the unit disk of the complex plane, and 
6q is arbitrary. Here, we have changed the sign with respect to the original Garrison- Wong 




(A.1) 
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paper and introduced arbitrary ^o- The inner product in is defined by 

eo+2TT 



(glf) = / d^/(e-'^)/(e-'^). (A.2) 

00 

The boundary value of / is given by a convergent Fourier series, 

-I oo 

fie-'') = ^Ecne-''^^ (A.3) 

which does not contain coefficients c„ with negative n. 

Subsequently, Popov and Yarunin [1973] established the connection of this operator 
to the Susskind and Glogower [1964] exponential phase operators E± of the form: 

= 0^ + + i [ln(l - e'^° E+) - ln(l - e"'^" E_)] . (A.4) 

The operators and £"+ (E-)^ are defined by the annihilation and creation operators 
d and of the mode: 

i_ = (a^a + l)-^/^a, i+ = (a^a + l)-^/^ [^_, ^+] = |0)(0|, (A.5) 

where |0) is the vacuum state [E_ is another notation for the Susskind-Glogower expo- 
nential operator (2.5)]. 

Let us consider the "phase states" 



1 oo 

1^) = j:^^PM)\n), (A.6) 



which are the right and left eigenstates of the operators E^ and E^: 

E.\e) = 

{e\E^ = eM-^0){e\. (A.7) 
The states (A.6) are not orthogonal but allow for the resolution of the identity operator: 

j de\e){e\ = i. (a.s) 

With the aid of eqs. (A.7) and (A. 8) to the operator, eq. (A.4) can be rewritten in the 
form (Bergou and Englert [1991]): 

= ^0 + TT + i j d9\e){e\ {ln[l - e-'(^-^°)] - ln[l - e'(^-^°)]} 

= j de\e)e{e\. (a.9) 



00 
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Since the states (A. 6) are not orthogonal, they are not eigenstates of the Garrison- Wong 
phase operator. 

From eq. (A. 9), we have 



(^I*GW|/) 



j de{g\e)e{e\f). 



(A.IO) 



Taking the field states |/) in the form 



I/) - E 



we then have 



Cn, n) 



n=0 



-j^ oo 



n=0 



(A.ll) 



(A.12) 



which has the same form as eq. (A. 3), and we can consider the phase operators (A.l) 
and (A.4) as equivalent. However, we should keep in mind that the Garrison- Wong phase 

operator is defined on a dense set of state vectors, which for mathematical consistency 
and the requirement that the number-phase commutator should be —i, imply /(—I) = 0. 
Unfortunately, when approximating even simple physical states on this dense set, one 
finds rather undesirable properties (Bergou and Englert [1991]). 
Since the states (A. 6) are not orthogonal, we have 



^GW 7^ / d0 0''\0){0\ {k>l), 



(A.13) 



and for the expectation values 



9o+2tt 

(/I^Gwl/) 7^ / d0 0'\{0\f)\' {k>l). 



(A.14) 



This means that the quantity |(6'|/)p cannot be interpreted as a phase distribution func- 
tion. To find the Garrison- Wong phase distribution function, we must calculate the quan- 
tity: 



Pgw{0) = \Gw{0\f)\' 



E Gw{0\n) 



n=0 



(A.15) 



where the vector |^)gw is the eigenvector of the Garrison- Wong phase operator. The 
function gw(^|'^) has a quite complex structure (Garrison and Wong [1970], Popov and 
Yarunin [1973, 1992]), but it can be found from the recursive formulas given by Garrison 
and Wong [1970], which are: 



G 



w{d\n) 



1 . e-00 

— sm — - — 

TT V 2 



.1/2 



MO), 



(A. 16) 
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where, for n > 1: 



MO) = -El-- ln-m{0)MO), 
m=0 ^ ^ 



(A. 17) 



ln{0) 



27r 



-2n 



d9' In 1^' - ^1 e'"^' - 



2n 



(A.18) 



-70 (6) 



(A.19) 



loid) - ~ + ^ [(27r + ^0 - ^) ln(27r + - ^) + - ^o) ln(^ - ^o)] ■ (A.20) 

The formulas (A.16)-(A.20) were used by Gantsog, Miranowicz and Tanas [1992] to cal- 
culate the Garrison- Wong phase distribution for some real states of the field showing that 
their symmetry is incompatible with the symmetry of the phase distributions obtained 
from the Pegg-Barnett as well as the s-parametrized phase approaches. In the Garrison- 
Wong approach even vacuum has a preferred phase, which is hardly acceptable on physical 
grounds. The recursive relation (A. 17) has the following solution (Miranowicz [1994]): 



MO) 



-70 (9) 



{ni,mi} t=l 



rrii 



-i::i&). 



(A.21) 



where the sum over {nj,mj} is taken under the condition Y^\=i'^i'<^i — ^ ^iid after inte- 
gration the functions 7^(6') [eq. (A.18)] take the form: 



in{e) 



27rin 



In 



27r 



ITT 



+ e'"^(Ei[in(27r + eo -^)] -Ei[-in(^-^o)])} (A.22) 

in terms of the exponential integral Ei(x). Equations (A.21)-(A.22) are more convenient 
for numerical calculations than eqs. (A.17)-(A.18). 

Substituting cq. (A. 6) into eq. (A. 9) and performing the integration over 9 yields the 
following number-states expansion for the Garrison- Wong phase operator [compare to 
eq. (2.14)]: 

exp[i(n — n')^o] \n){n'\ 



GW = 6*0 + TT ^ 



1 [n 



TV 



leading to the number phase commutator 

[$GW,ata] = -i(l-27r|^o)(^o|), 



(A.23) 



(A.24) 



and for the states for which (^o|/) = 0, the second term on the right-hand side vanishes, 
giving the value demanded by Garrison and Wong [1970]. A detailed comparison of the 

Garrison- Wong and Pegg-Barnett formalisms was given by Barnett and Pegg [1992] and 
by Gantsog, Miranowicz and Tanas [1992]. The difference between the two formalisms 
is, in mathematical sense, the difference between the weak and strong limits for the phase 
operators that is taken when cr — >■ oo (Vaccaro and Pegg [1993]). 
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Appendix B. States for the Pegg-Barnett Phase For- 
malism 



The Pegg-Barnett optical phase operator (2.12) is constructed in a finite (a + 1)- 
dimensional Hilbert space 'H'^'^^ spanned by the number states |0), |1), . . . , |cr). Hence, 
all other quantities, such as states, operators or probability distributions, analyzed within 
the Pegg-Barnett formalism, should also be defined in the same (cr-l- l)-dimensional state 
space l-i^"^^. Buzek, Wilson- Gordon, Knight and Lai [1992] emphasized that it is not 
strictly correct to apply the definition (2.23) of the finite-dimensional phase distribution 

P{em) = \ia){em\f)? (wrOUg), (B.l) 

for the infinite-dimensional state 

oo 

I/) ^ l/){oo) = E^nln). (B.2) 

n=0 

The problem of the precise definition of states in T^f'^^can be overcome by assuming that a 
is large enough so that the differences between the states in the finite, 1-L^'^\ and infinite- 
dimensional, "H, spaces can be arbitrarily small in the sense of the Cauchy condition (Pegg 
and Barnctt [1989]): 

V =1 " 

l-Y.\cn? < e. (B.3) 

^ ^ n=0 



The precise finite-dimensional phase 
Gordon, Knight and Lai [1992]): 

P{e^) = 

for the (cr -|- l)-dimensional state 
which is properly normalized, 

{a){f\f){a) -- 



distribution reads as follows (Buzek, Wilson- 

MOm\f)ia)\' (B.4) 



Ec^V), (B.5) 

n=0 



n=0 



for arbitrary a. The main problem resides in the construction of the normalized {a + 1)- 
dimensional states \f){a)- We restrict our attention to finite-dimensional coherent states 
only. However, other finite-dimensional states of the electromagnetic field can be defined 
in a similar manner; e.g., squeezed states (Buzek, Wilson- Gordon, Knight and Lai [1992]), 
even and odd coherent states (Zhu and Kuang [1994]), phase coherent states (Kuang and 
Chen [1994a, 1994b], Gangopadhyay [1994]) and displaced phase states (Gangopadhyay 
[1994]). 

There exist several generalizations of coherent states comprising the finite-dimensional 
case (see Zhang, Feng and Gilmore [1990] and references therein). It is possible to define 
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coherent states using the concept of Lie group representations (see, e.g., Pei^ina, Hradil and 
Jurco [1994]), or to postulate the validity of some properties of the infinite-dimensional 
Hilbert-space coherent states for the finite-dimensional coherent states. We present two 
definitions of the latter case. Firstly, the coherent states |Q;)(o^)in (a + l)-dimensional 
Hilbert space of a harmonic oscillator can be defined in the Glauber sense by the action of 
an analogue of the Glauber displacement operator D^'^\a) on the vacuum state |0) (Buzek, 
Wilson-Gordon, Knight and Lai [1992]): 

= D^''\a)\0) = exp(aa^ -Q;*a)|0). (B.7) 

The operator D^'^^ (a) is given in terms of the modified annihilation operator: 

d = exp(i<l>e)\/i^ 

= \Q){l\ + V2\l){2\ + --- + V^\a-l){a\ (B.8) 

and modified creation operator ci^ The coherent states |a)(o.)are close analogues of 
Glauber's (i.e., infinite-dimensional) coherent states \a). They were introduced and 
discussed by Buzek, Wilson- Gordon, Knight and Lai [1992], and their analytical Fock 
expansion was found by Miranowicz, Pi§,tek and Tanas [1994] in the form [eq. (B.5)] 
I/) — I a), with the superposition coefficients: 

Vn! T+1 k=o 

Here, xi = Xi^'^^^ are the roots of the modified Hermite polynomial of order (cr -|- 1), 
Re^+iixi) = 0, He„(x) = 2""/^H„(a;/y2), and a = |q;| exp(i^). 

Kuang, Wang and Zhou [1993, 1994] defined the normalized finite-dimensional co- 
herent states in another manner by truncating the Fock-basis expansion of the Glauber 
infinite-dimensional coherent states or, equivalently, by the action of the formally designed 
"displacement" operator exp(aa^) exp(— So) on the vacuum state. This approach is close 
to that of Vaccaro and Pegg [1990b] in the construction of a finite-dimensional Wigner 
function for coherent states. The states |a)((T) can be defined as follows (Kuang, Wang 
and Zhou [1993]): 

= A/'^'^^exp(aat)|0) = Ec^l^), (B-IO) 

n=0 

where 

c(r) = ^f('^^^. (B.ii) 

and the normalization constant is (Opatrny, Miranowicz and Bajer [1996]): 

^^"^ = {(-l)"L,^"i(|an}~'^' (B.12) 
in terms of generalized Laguerre polynomials L^(x). 
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The differences between the finite-dimensional coherent states (B.7) and (B.IO) were 
discussed in detail by Opatrny, Miranowicz and Bajer [1996] using the finite-dimensional 
Wigner function (Wootters [1987], Vaccaro and Pegg [1990b]) and in terms of the Stokes 
parameters. 

In the limit a — >■ oo, the coherent states |Q;)(o.)and go over into (a = a): 

lim \oi){a\ = lim \a)(cr-\ = la), (B.13) 

as was shown analytically by Opatrny, Miranowicz and Bajer [1996]. However, the states 
|a)(o-)and |a)(CT) are essentially different, particularly for |q;|^, |a|^ > ex, from the ordinary 
(infinite-dimensional) Glauber coherent states \a) as revealed by their photon- number, 
squeezing and phase properties (Buzek, Wilson-Gordon, Knight and Lai [1992], Kuang, 
Wang and Zhou [1993, 1994], Miranowicz, Pi§,tek and Tanas [1994]). Let us only mention 
that the well-known property of the ordinary coherent state jo;) for the mean photon- 
number is not fulfilled in the case of the finite-dimensional coherent states: 



(<,)(a|n|a)(^) 
(^) (a|77,|a)((^) 



^ {a\h\a) = \a\\ (B.14) 



The finite-dimensional states discussed here are not only mathematical structures. 
A framework for their physical interpretation is provided by cavity quantum electrody- 
namics and atomic physics. Moreover, they can be generated, e.g., in a single- mode 
resonator. Several methods have been proposed for the preparation of an arbitrary field 
state (e.g., Vogel, Akulin and Schleich [1993], Garraway, Sherman, Moya-Gessa, Knight 
and Kurizki [1994] and references therein), which can be readily applied to prepare these 
finite-dimensional states. Recently, Leohski and Tanas [1994] have presented a scheme 
of field generation in a cavity containing a nonlinear Kerr medium,kicked periodically 
with classical pulses. The field generated in this process is actually the finite-dimensional 
coherent state |Q!)(CT)in Hilbert space %'^"^ol arbitrary dimension. 
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